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Abstract 

We perform a general analysis of the dynamic structure of two classes of relativistic lagrangian field theories 
exhibiting static spherically symmetric non-topological soliton solutions. The analysis is concerned with (multi-) 
scalar fields and generalized gauge fields of compact semi-simple Lie groups. The lagrangian densities governing the 
dynamics of the (multi-) scalar fields are assumed to be general functions of the kinetic terms, whereas the gauge- 
invariant lagrangians are general functions of the field invariants. These functions are constrained by requirements of 
regularity, positivity of the energy and vanishing of the vacuum energy, defining what we call "admissible" models. In 
the scalar case we establish the general conditions which determine exhaustively the families of admissible lagrangian 
models supporting this kind of finite-energy solutions. We analyze some explicit examples of these different families, 
which are defined by the asymptotic and central behaviour of the fields of the corresponding particle-like solutions. 
From the variational analysis of the energy functional, we show that the admissibility constraints and the finiteness 
of the energy of the scalar solitons are necessary and sufficient conditions for their linear static stability against small 
charge-preserving perturbations. Furthermore we perform a general spectral analysis of the dynamic evolution of 
the small perturbations around the statically stable solitons, establishing their dynamic stability. Next, we consider 
the case of many-components scalar fields, showing that the resolution of the particle-like field problem in this case 
reduces to that of the one-component case. The study of these scalar models is a necessary step in the analysis of 
the gauge fields. In this latter case, we add the requirement of parity invariance to the admissibility constraints. 
We determine the general conditions defining the families of admissible gauge-invariant models exhibiting finite- 
energy electrostatic spherically symmetric solutions which, unlike the (multi-) scalar case, are not always stable. The 
variational analysis of the energy functional leads now to supplementary restrictions to be imposed on the lagrangian 
densities in order to ensure the linear stability of the solitons. We establish a correspondence between any admissible 
soliton-supporting (multi-) scalar model and a family of admissible generalized gauge models supporting finite-energy 
electrostatic point-like solutions. Conversely, for each admissible soliton-supporting gauge-invariant model there is an 
associated unique admissible (multi-) scalar model with soliton solutions. This shows the exhaustive character of the 
admissibility and stability conditions in determining the class of soliton-supporting generalized gauge models. The 
usual Born-Infeld electrodynamic theory and its non-abelian extensions are shown to be (very particular) examples 
of one of these families. 
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1. Introduction 

In the context of field theory the interest in extended configurations describing fields 
associated to point-like particles (particle-like solutions) dates back to the works of Mie [I] 
and Born and Infeld (B-I) [2], [3], aimed to solve the problem of the divergent self-energy of 
the electron field in Classical Electrodynamics. Today, this kind of problems fall inside the 
large domain of Soliton Physics, whose methods and applications concern most branches of 
physical sciences. 

Although there is not an universally accepted definition of the concept of soliton, we shall 
adopt here minimal defining properties which seem to be widely assumed in most contexts: 
Solitons are stable, finite- energy solutions of conservative non-linear differential equations. 
However, the accepted meaning and content of the term stable is not universal. In a strong 
sense, it refers to the existence of soliton entities which can be identified (if present) in field 
configurations and are preserved by the dynamic evolution of the system. With this defini- 
tion, the analysis of such configurations in terms of many solitons, interacting via radiative 
field exchanges, becomes possible. This kind of stability arises in some field theoretical mod- 
els (most of them in one-space dimension) exhibiting topological conservation laws, related 
to a non-trivial structure of the vacuum |3] . In these cases the conserved topological charges 
identify the presence of the topological solitons and ensure their preservation. Examples of 
topological solitons in three-space dimensions are the monopole of 't Hooft and Polyakov [5] 
or the Skyrmion [6j. In a weak sense, stability is identified with linear stability, i.e. with the 
preservation of the soliton identity against a certain class of small perturbations for which 
the soliton configuration is a minimum of the functional of energy. This restricted class 
of perturbations is defined through boundary conditions which amount, in general, to the 
preservation of (non-topological) charges associated with the soliton. The conservation laws 
of these charges may be implicitly contained within the structure of the field equations (as 
in electrodynamics) or be consequences of constraints imposed on the external sources, to 
which the field is coupled. Here we shall be only concerned with this kind of non-topological, 
finite-energy, weakly-stable soliton solutions of local relativistic lagrangian field theories in 
three-space dimensions. 

Non-linear field theories supporting soliton solutions have been widely studied and applied 
in several contexts of Theoretical Physics. Let us give some examples: 

1) In the context of elementary particle physics, let us mention the phenomenological 
description of the hadron structure and the hadronic interactions in terms of topological 
solitons in the framework of the Skyrme model [6j. Other descriptions of this structure in 
terms of non-topological solitons have been performed through the effective approach to the 
low-energy regime of Quantum Chromodynamics (QCD) given by the Friedberg-Lee model 
[7], [8] and related theories. Let us also mention the generalized chiral-invariant model of 
Deser, Duff and Isham [9j, whose lagrangian density is a rational power (3/2) of the non- 
linear sigma model lagrangian, chosen in order to circumvent Derrick's theorem [IQl. This 
model and its extensions support topological solitons [11]. In the same context let us also 
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mention the toroidal solitons of Ref . [12] , which might describe glueball collective states in 
the low-energy limit of QCD, as suggested in Ref. [T3] . 

2) In nuclear physics, the analysis of high-density hadronic matter and its chiral phase 
structure has been performed in terms of systems of skyrmions in the framework of effective 
field theories implemented with the large A^^^ behaviour of QCD |14j . 

3) The glueball collective states, as soliton solutions of non-abelian B-I gauge field models 
[T5] . The introduction of these generalized gauge models supporting soliton solutions and 
their extensions to higher dimensions is suggested by string theory, since some of them arise 
in the low-energy limit of D-Branes [16] . Moreover, it has been shown that the B-I extension 
of the basic lagrangian of the Skyrme model leads to stable soliton solutions. This extension 
was proposed as an alternative to the inclusion of the quartic "Skyrme term" in the original 
lagrangian of the non-linear sigma model [T7]. This "ad hoc" term stabilizes the soliton by 
preventing Derrick's scaling, but its introduction is not justified from more fundamental 
reasons. 

4) In the last two decades there has been an increasingly amount of works on self- 
gravitating field configurations. The aforementioned presence of B-I actions in the low-energy 
physics of D-branes, whose fundamental excitation is gravity, is one of the motivations for 
this renewed interest (see e.g. [H]). But the search for self-gravitating field configurations, 
as solutions of the Einstein equations for gravity coupled to different kinds of fields, is an 
older topic [Hj. Indeed, in a four- dimensional flat space-time. Derrick's theorem [TO] and 
several non-existence theorems [20] limit drastically the class of lagrangian field theories 
supporting soliton solutions. However, the coupling to gravity can remove these obstruc- 
tions and thus allow for particle-like solutions in some cases. One example of this is the 
pure Yang-Mills theory, which does not support glueball solutions in Minkowski space [20] 
but exhibits particle-like solutions in curved space-time [21]. Another example is the regu- 
lar black- hole solution of Einstein's equations coupled to a non-linear electrodynamics [22] . 
Let us also mention that theories supporting soliton solutions in fiat space, as the abelian 
and non-abelian Born-Infeld models, have been extended to curved space leading to black- 
hole-like soliton solutions [23]. Finally, we cite the soliton stars introduced by T.D. Lee as 
self-gravitating coherent quantum states with the features of non-topological solitons [5].[2^. 

5) In the fast-evolving context of modern Cosmology let us mention some problems for 
which non-linear field theories supporting soliton solutions have been invoked. It has been 
suggested that time-dependent but non-dispersive scalar solitons (Q-balls [25]) may account 
for the behaviour of self-interacting dark matter [26] . There is also the suggestion that non- 
topological solitons might have been formed in a second-order phase transition in the early 
Universe, and contribute significantly to the present mass density [27] . As another example, 
the Born-Infeld generalization of SU(2) non-abelian gauge field theory, coupled to tensor- 
scalar gravitation, has been used for the description of dark energy [2S]. Also scalar field 
models with lagrangian densities which are generalized functions of the kinetic term have 
been used to drive inflationary evolution in the early Universe (k-infiation) [2^]. Solitonic 
configurations of these k-essence fields have been used to reproduce some properties of dark 
matter as well [30] . 

All these considerations underline the interest of non-linear lagrangian field theories and 
their eventual soliton solutions, mainly for the cases of (one and many-components) scalar 
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field models and generalized gauge- invariant field theories. By generalized gauge field theories 
we mean models for gauge fields of compact semi-simple Lie groups, with lagrangian densities 
defined as general functions f{X, Y) of the two standard first-order gauge invariants, namely 
X = tr{Ffj,y- F^^) and Y = tr{F^y- F*^^). Aside from the already mentioned Born-Infeld-like 
gauge models, defined by the very particular B-I choice of this function, there are not in 
the literature systematic studies on solitons for general gauge-invariant field theories. It is 
one of the main purposes of this work to perform a study on this class of theories in three 
space dimensions, by setting the conditions on the lagrangian functions determining the 
families of these models which are physically admissible and support elementary solutions 
which are non-topological solitons lJ. The requirements for admissibility adopted here refer 
to the positive-definite character of the energy, the vanishing of the vacuum energy and 
the regularity, uniqueness and definiteness in all space of the elementary solutions. These 
conditions are targeted to deal with physically meaningful theories. 

In principle, generalized gauge-invariant lagrangians are candidates to describe the dy- 
namics of the gauge fields in gauge theories of fundamental interactions. If one accepts 
the fundamental character of string theory and the aforementioned results, referenced in 
[16], the description of the gauge-field dynamics through some generalized lagrangians, re- 
garded as effective field models of string theory, could be more "fundamental" than the usual 
Maxwell-like choice ip{X, F) ~ X, a "minimal prescription" which should be understood as 
a low-energy (or weak- field) approximation limit. Nevertheless, from the field-theory point 
of view, this minimal prescription is generally assumed to describe the fundamental dynam- 
ics of the gauge fields in the gauge-invariant lagrangians with coupling to other (generally 
fermionic) sectors. In this case, when the high-energy degrees of freedom are integrated 
out in the path integral of the original action, generalized gauge-invariant models emerge 
as effective lagrangians, containing new non-linear self-couplings of the gauge fields which 
account, at a classical phenomenological level, for quantum effects and interactions with 
the removed heavy- mode sector [31]. Historically, the first example of this kind of effective 
lagrangians was obtained by Heisenberg and Euler [32] in the context of Quantum Electro- 
dynamics (QED). It accounts for the non-linear effects of the Dirac vacuum on low-energy 
electromagnetic wave propagation, calculated to lowest order in the fine structure constant. 
When higher order operators are included we are lead to a sequence of effective lagrangians 
which take the form of polynomials in the field invariants, arranged as an expansion in 
operators of increasing dimensions [33]. An interesting question arises here, related to the 
possibility that these effective lagrangians could exhibit soliton solutions, even though the 
bare lagrangian does not. In the following sections we shall give explicit examples for which 
the soliton elementary solutions of an effective model may be interpreted as finite-energy 
fields of point-charges screened by the vacuum effects, whereas the elementary field of the 
bare theory is energy-divergent. 

For scalar fields. Derrick's theorem [TU] imposes severe restrictions on the lagrangian mod- 
els supporting time-independent soliton solutions in three space dimensions. One of the ways 
to circumvent the hypothesis of Derrick's theorem is based on the choice of the lagrangian 
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density as a function of the kinetic term alone [31] (see Eq. fl4.ip below). This choice seems 
rather arbitrary from a physical point of view. Nevertheless, it is the natural restriction for 
scalar fields of the lagrangian densities of the generalized gauge-invariant models. Moreover, 
as we shall see in the sequel, many of the necessary results in the analysis of the soliton 
problem in generalized (abelian and non-abelian) gauge-invariant theories (characterization 
of the families of soliton-supporting lagrangians, explicit determination of such models, anal- 
ysis of stability etc.) can be obtained from similar results more easily established for scalar 
models. Consequently, the detailed study of these scalar models will take an important place 
in this work. 
The paper is organized as follows: 

In section 2 we consider the scalar field models in detail. Many general results concerning 
the properties of these families of scalar field theories and their associated soliton solutions 
can be obtained without the specification of the explicit form of the lagrangian functions. 
After an initial discussion on the admissibility conditions to be imposed on these functions, 
we solve the field equations for SSS solutions of generic models. This is achieved by obtaining 
the generic expression of a first-integral, which allows the determination of the field strength 
once the form of the lagrangian function is specified. We analyze the expression of the 
integral of energy for these solutions and determine the conditions that must be satisfied 
(at r = and as r ^ oo) in order for them to be of finite-energy. These conditions imply 
supplementary restrictions to be imposed on the lagrangian functions which allow for an 
exhaustive characterization of the admissible models supporting finite-energy SSS solutions. 
We shall leave for section 6 the study of the stability of the elementary solutions for these 
and other models. 

In section 3 we introduce explicit examples of the different families of admissible scalar 
models supporting SSS soliton solutions. The first example is a large family of polynomial la- 
grangian functions including the scalar versions of the Euler-Heisenberg effective lagrangian 
of Electrodynamics and the sequence of higher-order effective corrections. From these ex- 
amples and the results of section 5 on generalized electromagnetic field theories it can be 
shown that the sequence of effective lagrangians for Electrodynamics describing low-energy 
photon-photon interaction, obtained in the perturbative expansion [33] , exhibit finite-energy 
point-like solutions. The second example is a two-parameter family generalizing the scalar 
version of the Born-Infeld model. Both the first and second examples exhibit soliton solu- 
tions which are asymptotically coulombian and differ by their behaviour at the center. The 
third example is a three-parameter family of models behaving like the B-I one at the cen- 
ter, but exhibiting diverse asymptotic behaviours. The fourth example is a three-parameter 
family that supports exponentially-damped soliton solutions. When properly extended to 
the generalized gauge-invariant case we obtain examples of gauge theories where non-linear 
self-couplings lead to short-ranged interactions without any symmetry breaking mechanism. 

In section 4 we extend the results of the analysis of one-component scalar fields to the case 
of A^-components scalar fields. For a given form of the lagrangian density as a function of 
the rotationally-invariant kinetic term (X = J2iLi d^(j)id^4>i) the N components of the SSS 
solutions have the same form, as functions of r, as the SSS solution of the one-component 
model associated to the same form of the lagrangian density function. The scalar charge of 
the associated one-component solution corresponds to the mean-square of the N external 
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scalar charges associated to the A^-component solutions and the energy (when finite) is the 
same in both cases. There is thus a degeneracy of the energy on spheres of the charge space 
in the A^-component case. These results will be useful in the subsequent analysis of the 
generalized non-abelian gauge-invariant models. 

In section 5 we consider the abelian and non-abelian generalized gauge-invariant theo- 
ries. After defining the dynamic problem and the admissibility conditions, we analyze the 
equations for ESS fields. We prove that the solutions of these equations, in the abelian 
and non-abelian cases, can be built from SSS solutions of the one and many- components 
scalar field problems, respectively. As a consequence there is a correspondence between scalar 
models and families of generalized gauge models in such a way that the SSS solutions and 
the corresponding ESS solutions have the same functional form. Moreover, if the energy 
of a SSS solution is finite, so is the energy of the ESS solutions of the corresponding gen- 
eralized gauge family. The results of section 2 characterizing the admissible scalar models 
with finite-energy SSS soliton solutions, characterize also the admissible generahzed gauge 
models with finite-energy ESS solutions, but the latter ones are not always stable. Stability 
requires now supplementary conditions to be satisfied by the lagrangian densities, which will 
be determined in section 6. 

Section 6 is devoted to a detailed analysis on the linear stability of the elementary solu- 
tions of scalar and gauge models. This analysis leads to necessary and sufficient conditions 
for the stability of the soliton solutions, going beyond the necessary conditions demanded 
by Derrick's theorem. For (multi-) scalar models the variational study of the functional of 
energy and the spectral analysis of the small perturbations around the SSS solutions prove 
the static and dynamic stability of all finite- energy SSS solutions of admissible lagrangians. 
For admissible abelian and non-abelian generalized gauge models, a similar analysis allows 
to determine supplementary (necessary and sufficient) conditions to be satisfied by the la- 
grangians in order for their finite-energy ESS solutions to reach stability. 

We conclude in section 7 by drawing some perspectives and future developments. 

Our analysis deals with fields in three-space dimensions, but most of our results can be 
straightforwardly generalized to other spatial dimensions. 

Short reports of some of the main results developed here have been already published in 



2. Scalar solitons 

We begin with lagrangian densities for scalar field potentials 0(x), defined in a four di- 
mensional Minkowski space-time as 



where f{X) is a given continuous and derivable function in the domain of definition {Q) 
which is assumed to be open, connected and including the vacuum {X = d^cf) ■ d^cf) = 0). 
For future purposes, we also require f{X) to be monotonically increasing (more precisely, 
df/dX > 0,VX 7^ e n and df/dX > for X = 0) and df/dX to be continuous for 
X < 0{X E Q). In absence of a coupling to external sources, the associated field equations 
take the form of a local conservation law 



Refs. ^ and 



L = f{d,4>-d^<P) 




6 



d^Jf" = 0, (2.2) 
where the conserved current is 

=f iX)dy, (2.3) 

with / (X) = df/dX. In these models, the D'Alembert hnear wave equation corresponds to 
f{X) = X/2. For the SSS solutions 0(r), Eq.([22D has the first-integral 

/ (-0'') = A, (2.4) 

where 0' = dcp/dr, and A is the integration constant. This is an algebraic equation which 
allows, in principle, the determination of the field strength as a function of r and A. The 

positivity of / {X) requires both A and (f)'{r) to be simultaneously either positive or negative. 
We can then consider only the positive-sign case without loss of generality. Strictly speaking 
Eq. (12.41) determines the field (p (r) only for r > 0. If we replace the solutions of Eq. (l2.4p 
in Eq. (12.21) we do not obtain zero, but a Dirac 6 distribution of weight 47rA. We can then 
identify this parameter with the central scalar charge source of the (at rest) SSS solution, in 
analogy with the point-like charges in the Maxwell theory. Alternatively, in some cases as, 
for instance, the non-linear electromagnetism of Born-Infeld, this charge may be interpreted 
as a continuous charge density distribution in space. For non-linear electromagnetic models 
the continuous interpretation of the charge is rather natural, owing to the conservation of the 
electric charge as a consequence of the field equations, but this is not so for the scalar models 
[37] . Nevertheless, following the analogy with the electromagnetic case, we can define for the 
models (12. ip the total scalar charge associated with a given static asymptotically vanishing 
field solution (f){r) as 



f (^)V(0) 



I f (X)V(0) ■ da, (2.5) 

5" 



which, owing to the field equation (12. 2p . vanishes for everywhere-regular solutions. For the 

SSS solutions of (12.40 we have / (X)V(0) = A^ and the total charge equals A. We can 

then define the spatial charge-density distribution as cr(r) = (l/47r) / (O)V^0, which gives 
for the total scalar charge of the SSS solutions 

^ f dsff (0)VV = Jhn / (O)rV(r) = Jim (2.6) 

• 

Clearly, this interpretation is only possible if / (0) is finite, or equivalently, if the function 
(r) goes to a constant as r — ^> 00 (asymptotically coulombian fields). This function must 
also vanish as r — 0. As we shall see at once, this latter condition is fulfilled for all models 
with finite-energy SSS solutions, but the former defines a sub-class of those models (see case 
B-2 below in this section) to which the scalar version of the Born-Infeld model belongs. 



^ As in this linear case, plane waves of the form (j> = (j>(k^-x^^), with fe^ = k^k^^ = 0, are solutions of Ea.l l2.2l l but superpositions 
of such waves are not, in general. 
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Once the form of f{X) is fixed, equation ( ]2.4p gives the expression of the field 0'(r) = 
(r, A) in imphcit form and allows the determination of the potential 0(r) (up to an additive 
arbitrary constant) through a quadrature. We also note that Eq.f l2.4p implies that (r), if 
single-branched E], is necessarily a monotonic function of r. Moreover, the solution depends 
on r and A through the ratio r / \fX. This is a straightforward consequence of the invariance 
of the solutions of the field equations (12. 2p under space-time scale transformations. Indeed, 
if 0(?^, t) is a solution of this equation, so is the modified function 

(^(f,t,A) = A"V(Ar,At), (2.7) 

A being a positive constant. This is a symmetry of the solutions of the field equations without 
sources, but not an invariance of the action [38] . 
The canonical energy-momentum tensor associated to the lagrangian (12.11) is 

T^, = 2 / (X)5^09,0 - /(X)r/^„ (2.8) 
and the energy density in the SSS case becomes 

p = T°o = -/(-0'2), (2.9) 
whereas the total energy is 

oo 

e(A) = -47r / rV(-0''(r, k))dr = A^/h{A = 1), (2.10) 







the last equality being a consequence of the aforementioned scale invariance. Using the first- 
integral (12.41) and integrating by parts we obtain the following useful expression for the 
energy in this SSS case 



e(A) = ^|[0(r,A)]| 



oo 



r0'(r,A) + ^/(-0'2(r,A)) \. (2.11) 




As we shall see, if the energy of the SSS solutions is finite the second bracket in the r.h.s. 
vanishes and this expression reduces to 

6(A) = ^[0(oo,A)-0(O,A)], (2.12) 

which shows that the potential 0(r) must be a bounded function of r, defined up to an 
arbitrary constant. Note also that Eq. (l2.12l) has the form of the potential energy of a point- 
like scalar charge of value 2A placed at infinity in the soliton field. 

There is another way round to obtain the expression (12.12^ . Indeed, by performing the 
usual reescaling of Derrick's theorem [10], i.e. by defining a uniparametric family of trans- 
formations 

0a(x) = (f){Xx), (2.13) 



^ In some cases Eg. 112.41 1 can lead to discontinuities or several branches for the function ij> (r). We shall regard such cases 
as "unphysical" and rule them out from this analysis, considering only models for which the fields of the SSS solutions are 
(for r > 0) continuous, single-branched functions defined everywhere. We shall establish at the end of this section that the 
corresponding admissibility condition for the lagrangian densities is the strict monotonicity. 
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it can be easily checked that the condition of extremum of the energy against these reescahngs 
^M^^xYi |^^_^_ leads automatically to the relation fl2.12p when the condition of finiteness 
of the energy is assumed. Although Derrick's theorem is only a necessary condition for sta- 
bility, this signals the presence of a connection between stability and finite-energy condition 
of the SSS solutions considered here. This connection will be precisely established in section 
6. 

Going beyond the results of reference [31], we shall determine general conditions to be 
imposed on the functions f{X) in order to obtain physically consistent field theories, whose 
associated SSS solutions be stable and their energy (12.101) be finite (non-topological soli- 
tons). We first summarize some criteria of physical consistency adopted for the purposes of 
the present study (defining what we shall call "admissible" field theories) and obtain the 
associated restrictions on the lagrangian densities. Next we shall obtain the conditions for 
such admissible models to support SSS soliton solutions. 



2.1. Conditions on the energy functional 



If any of these models are to be used for descriptions of quantum physical systems, the 
possibility of their quantization becomes important. This implies supplementary conditions 
to be satisfied by the lagrangian densities, as necessary for any quantum extension. Thus, 
aside from the above mentioned continuity and derivability conditions, we must require the 
function f{X) to be defined everywhere (f2 = 3?), in order to allow the proper definition of 
the associated path integral. Here we shall call such models class- 1 field theories 0- 

Alternatively, as discussed in the introduction, we can consider these scalar models (and 
their generalized versions proposed below) as effective classical lagrangians of more "funda- 
mental" theories, including integrated high-energy and quantum effects through new non- 
linear couplings. In these cases we can relax the everywhere definiteness conditions of the 
admissible lagrangians and require their regularity only within a restricted domain of defini- 
tion (f2 C 3?), which is assumed to be open, connected and including the vacuum (0 G fi). We 
shall call these models class-2 field theories. Using this criterion, models such as the Born- 
Infeld one become admissible field theories belonging to this class. Although these models are 
essentially classical, quantum corrections to their particle-like solutions can be obtained by 
quantizing the field of small fluctuations around these ground states. Such fluctuations obey 
Euler-Lagrange linear field equations of admissible lagrangian densities defined everywhere 
(see Ref. [39] and section 6, Eqs. (l6.1ip and (I6.12p below). 

Obviously, a second condition to be imposed for admissibility in all cases concerns the 
positive definite character of the energy, which is required to hold in the entire domain of 
definition (fi) of the lagrangian density. The expression for the energy density in terms of 
f{X) is 



* These conditions exclude models such as the scalar version of Born-Infeld Electrodynamics [S]- In such cases it is always pos- 
sible to generalize the model, by continuing the lagrangian density function to the undefined regions through some prescription 
which must preserve the classical dynamical content of the initial model. Then, the quantum behaviour of the extended model 
would depend on this prescription. But this procedure will necessarily enlarge the space of solutions of the classical theory. 
As we shall show below, these extensions introduce new branches for the SSS solutions which become spurious at the classical 
level. We shall exclude such extended models from the present analysis since they are non-admissible according to our physical 
criteria. 
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p = r°o = 2 / (X) (^^^ - /(X) = 2X / (X) - /(X) + 2 / (X) (V0)' . (2.14) 
For the D'Alembert lagrangian we have /(X) = X/2, and the energy density reduces to 



P 



dt 



V0 



(2.15) 



which is positive for any non-constant function 0(t, r) and vanishes in vacuum. In order to 
obtain a similar behaviour in the general case (12.141) (requiring also the energy density to 
vanish in vacuum) we are lead to the necessary conditions 



/(0)=0 
/ (X) > vx. 



(2.16) 



For minimal sufficient conditions on /(X) let us analyze separately the cases X < and 
X > 0. On the one hand for X < the term (9^0)^ may vanish and the positivity of the 
energy density (12.14p requires 



/(X)<0 (VX<0). 



(2.17) 



On the other hand, for X > the term (V</))^ may vanish and the positivity of the energy 
requires 



p(X) > 2X / (X) - /(X) = X/(X) 



dX 



In 



> 0. 



f2.18) 



This equation, together with the conditions /(O) = and / (X) > (VX), lead to 

> 0, (2.19) 



in [l^ 



d 

dX ^ 

or, equivalently, the function (and hence /(X) itself) must be a positive monotonically 
increasing function of X for X > 0. Equation (12.19^ . together with the boundary condition 
/(O) = 0, fix the behaviour of /(X) around X = as 



/(X ^ 0+) ~ X^+", 
with a > —1/2, and the energy density behaves there as [3 
p{X 0+) > (1 + 2a)Xi+". 



(2.20) 



(2.21) 



° For values of a in Eqs. |(2?20] | and |(Z2T|| which he in the interval -1 < a < -1/2 the condition /(O) = is fulfilled, but 
• 

/ {X) diverges at X = in such a way that the energy density becomes necessarily negative in the neighbourhood of X = 0. 
For a < —1, f{X) diverges in vacuum. The limit case a = —1/2 is singular. 
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2.2. Conditions for finite- energy SSS solutions 



The convergence of the integral of energy f l2.10p for the SSS solutions is governed by the 
behaviour of the integrand near the limits r ^ oo and r — 0. This imposes supplementary 
conditions on the form of the function f{X) around the values of X{r) = — 0^(r) in these 
limits. Let us assume a power law expression for the field around these regions {4>' (r) ~ r'^ 
as r — i> oo or as r ~ 0); from the first-integral (12 ■4p we obtain the relation (valid for q ^ 0; 
for g = at r = see the discussion of the case A-2 below) 



df /dr 



2A0"(r) 



-2Agr«"^ 



(2.22) 



and in the limits of integration, /(r) behaves as 



fir) 



2Aq 



2-q 
if g 7^ 2, or as 



f{r) 4Aln(r) + D, 



(2.23) 



(2.24) 



if g = 2. The integration constants D in these expressions are easily related to the values of 
X and f{X) on the limits of the integral of energy Eq. (l2.10p . Around each of these limits 
the contributions to this integral take the form 



~ An J dr 

for q 2 and 
— An I dr 



2Ag „ o 
^-r« + Dr^ 



4Ar^ln(r) + Dr' 



(2.25) 



(2.26) 



for q = 2. Let us analyze separately the convergence of the energy integral around r ~ 
(case A) and in the asymptotic limit r — >■ cxd (case B). In case A the convergence of fl2.25p 
requires g > — 1 and we can distinguish three sub-cases: 

• A-1 ) If — 1 < g < the field (p (r) diverges at r ^ but the integral of energy converges 
there and the potential 0(r) is finite at the origin. Then, as r approaches zero, X —>■ —oo 

and f{X) and / (X) diverge as 



fix) 



q-2 



-X)^ 



f (X) ~ (-x)-^^ 



(2.27) 



(see figure 1). Such solutions can be stable and finite-energy SSS fields (depending on their 
large-r behaviour) and, in this sense, they might be considered as genuine non-topological 
solitons. 



^ Although this assumption excludes some transcendent behaviours such as the asymptotic exponential damping, our conclu- 
sions will remain valid for models exhibiting these "short- ranged" SSS solutions. In fact such models are included in the case 
B-3 below. Note that damped oscillatory behaviour at infinity is excluded by the monotonicity of the SSS field solutions. 
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B-3 




Fig. 1. Different possible behaviours of the admissible lagrangians supporting finite-energy SSS solutions. Note the existence 
of two A-2 cases, both ending at a finite value X = —C^ (corresponding to the maximum value of the field strength) but with 
bounded or unbounded values of f{X) there. A-1 branch is parabolic with divergent slope as X ^ — oo. 

• A-2 ) If g = the field (j)' (r) goes to a constant value at the origin (0'(O) = C, which 
corresponds to X + = 0) and can be written around this point as 

0'(r)~C-^r", (2.28) 

where C, 6 and a are positive constants [3. Then / (X) diver ges there as 

/ (X) ~ (X + C^)-^. (2.29) 
Consequently, for cr 7^ 2, /(X) behaves around X = —C^ as 

/W- ^ /| A +A, (2.30) 

where A is a constant. For a = 2 this behaviour becomes 

/(X) ~ln(X + C2). (2.31) 

We see that when a < 2 these lagrangians diverge at X = — (figure 1). However they 
can be accepted as admissible class-2 field theories if this point is located on the frontier of 
the (open) domain of definition where the lagrangian density must be regular (see the third 
example of the following section). Consequently the set X < —C^ must be excluded from 
the domain of definition (Q) in this case. This requires a to be an irrational number or four 
times the irreducible ratio between any natural and an odd natural. When a > 2 these values 
also exclude the same region, leading again to admissible class-2 field theories. However, in 

this case the lagrangians are finite in X = — (even if / (X) diverges there). Then, for 
rational values of a > 2 which are irreducible ratios of an odd natural and any natural 
0, the lagrangians are defined and continuous for any X around X = — C^, exhibiting a 



Note that the scalar version of the Born-Infeld model is an example which belongs to this case, corresponding to tr = 4. 
* Other values of a, for which the lagrangian is also defined for X < —C^, lead to negative energy densities there and 
consequently must be excluded. 
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vertical-slope inflexion in this point. If appropriately extended to all X, they lead to class- 
1 fleld theories which satisfy the everywhere positive deflniteness condition of the energy. 

Nevertheless they violate the requirement of continuity of / (X) for X < 0(g and this 
leads to associated point-like solutions exhibiting several branches. Indeed, as we shall see at 
once, this requirement for admissibility is introduced because it endorses the single-branched 
character of the SSS solutions. Consequently, in this case A-2 all admissible lagrangians must 
remain undeflned for X < — and thus belong to class-2 fleld theories. This implies (as 
happens in the Born-Infeld model [3]) the existence of a maximum value of the fleld strength 
(0'(r)<C). 

As a function of r, the energy density behaves around the center as 

^ ' 1-2/a ' ^ ^ 

for a ^ 2 and as 

rV(r) ~ 2r^ ln(r) + ln(2C^), (2.33) 
for a = 2. As expected, the energy integral converges there in both cases. 

• A-3 ) The case g > must be discarded. Indeed, in this case / (X) behaves as 

• 9+2 

/(X)~X" — , (2.34) 

around X = 0. Consequently, /(X) is singular in vacuum for < g < 2. For q > 2 the 
energy density for X — > 0"^ behaves as 

p(X) ~ --X^, (2.35) 

q 

and becomes negative around the vacuum (see also Eqs. (12.201) . fl2.2ip and the footnote there). 

In case B the convergence of f l2.25p in the r — oo limit requires q < —1 (in this case 
(r ^ oo) = and the integration constant D in (12.251) vanishes). Then the behaviour of 
/(X) around X = must beLj 

/(X)~X^. (2.36) 

Now the existence of the lagrangian on both sides around X = becomes crucial for the 
consistency of the theory and this imposes supplementary restrictions on the possible values 
of the parameter p. Indeed, the exponent in Eq. (12.361) must be the ratio of two odd naturals 
F°1 and this restricts the possible values of p to a sub-class of the rational numbers. Let us 
analyze separately three possibilities: 

• B-1 ) Consider flrst the case 1 < p < 2. We deflne P and Q as two positive odd natural 
numbers in such a way that P < Q and the ratio S = P/Q he irreducible. Then the 
admissible values of the exponent in (I2.36P are given by 

p + 2 3 



2p 2 + S ' 



(2.37) 



® For the sake of clarity we use here the parameter p = —q in the exponent, in terms of which <f>(r —> oo) ~ ^ with p > 1. 

If this exponent is the irreducible ratio between an even and an odd natural numbers the lagrangian is well defined on both 
sides of X = 0, but the energy density becomes negative for X < 0. 
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and the corresponding admissible values of p can be written as 

4 + 2E , , 

P=^3^- (2.38) 

Now / (0) = and the slope of the lagrangian vanishes in vacuum (see figure 1). 

• B-2 ) For p = 2 the lagrangian behaves around X = as the D'Alembert lagrangian 

f{X ^ 0±) ~ X, (2.39) 

and the soliton field becomes asymptotically coulombian (see figure 1). 

• B-3 ) For p > 2 the behaviour of the lagrangian is also given by Eg. (12.361) but now the 
admissible values of the exponent are constrained by 

^ = ^ > 1, (2.40) 

2p 1 + S 2' ^ ^ 

where S = P/Q must be the irreducible ratio between an even natural P and an odd natural 
Q such that Q > P. The corresponding admissible values of p are 

1 + S , , 

p = 2^-^. (2.41) 

As easily seen, the slope of the lagrangian diverges at X = in this case (see figure 1), but 
the energy density remains positive definite there. 

We conclude that the set of admissible models exhibiting finite-energy SSS solutions can 
be classified into six families which are the combinations of the cases A-1 or A-2, governing 
the central field behaviour and the cases B-1, B-2 or B-3, determining the asymptotic field 
behaviour. Moreover, any given scalar, monotonically decreasing, SSS function 0'(r), which 
satisfies boundary conditions of A-type at the center and of B-type asymptotically, is a finite- 
energy SSS solution of a particular lagrangian model belonging to one of these families. The 
explicit form of the corresponding lagrangian density can be found by integrating Eq. fl2.4p 
with respect to the variable X = —(f) ^(r) with the corresponding boundary conditions (see 
the fourth example of the next section). 



2.3. Conditions for stability 



We summarize here the main steps in the analysis of linear stability of the scalar SSS 
soliton solutions (the detailed calculations are given in section 6). The linear stability of these 
solutions requires their energy fl2.10p to be a local minimum against small charge-preserving 
perturbations. We consider finite-energy SSS solutions 0(r) and small static perturbations 
50(r), finite and regular everywhere and vanishing (as well as their radial derivatives) as 
r — > oo, in such a way that the scalar charge of the perturbed fields remains unchanged at 
the first order in the perturbations. For the static solutions of the field equations (12. 2p the 
first variation of the energy (I2.14p vanishes, while the second variation is positive if and only 
if the condition 

/ (X) + 2X 7 (X) > 0, (2.42) 

See Eq. 1 16. 7t and the footnote there. 
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is satisfied in all the range of values of X covered by the solution (X = — 0'^(r),O < r < 
oo). As we shall see this requirement is always fulfilled by the finite-energy SSS solutions 
of the admissible models defined in this section, proving their linear stability. A detailed 
spectral analysis of the small oscillations around these SSS finite-energy solutions has also 
been performed. It leads, for admissible models, to discrete spectra of eigenvalues and nor- 
malizable orthogonal eigenfunctions in Hilbert spaces, whose scalar products are built as 
three-dimensional integrals of their products with the functions f{X{r)) as kernels. In their 
temporal evolutions the perturbations to the soliton solutions remain bounded in this norm, 
confirming the stability (see subsections 6-1 and 6-2). 

2.4. Conditions for uniqueness of the elementary solutions 

We return now to Eg. (12.41) . which defines the SSS field solutions (j) {r), and analyze the 
conditions under which they are single-branched and defined everywhere functions. As al- 
ready mentioned, owing to the positivity of / {X) both A and (j) {r) must be simultaneously 
either positive or negative and then we can analyze only the case (r) > without loss of 
generality. Let us write (12. 4p under the form 

z{y) =yf (V) = ^ > 0, (2.43) 

where we have introduced the variable y = (f)' (r) > 0. The field strength function is given 

by the values of y obtained by cutting the curve z{y) with horizontal lines corresponding to 

the different values of r. Then, for the field (r) to be defined in all space, z must range 

from to oo and the uniqueness of the solution requires a single cut point on every z = 

constant line. This restricts z{y) to be a continuous monotonic function. As a consequence, 

the requirements of continuity (for X < 0) and strict positivity (for any X 7^ 0) imposed on 
• 

/ (X) at the beginning of this section, when establishing the admissibility conditions, are 
mandatory for the proper definition of the SSS solutions. Indeed, a glance at the form of 

the function z{y) reveals that the existence of a jump in the function / (X) at a finite value 
X < would lead to SSS solutions which are double-branched or undefined in some range 

of values of r. Moreover, if / (X) vanishes for a value Xq < (horizontal-slope inflexion 
point for /(X)), then z{y) vanishes at y = (— Xq)^/^ reaching a minimum there and the SSS 

solutions become necessarily multiple-valued (the vanishing of / (X) for a positive value 
Xq > is discarded by the energy-positivity condition (12.191) ). 

For soliton solutions the finiteness of the energy requires the origin to be a point of the 
curve z{y) (see figure 2) and the large- 2; behaviour is determined by the behaviour of the 
field at the center of the soliton. Then, in the cases where the field strength diverges as 
r — >• (case A-1 above) the uniqueness of the solution requires the curve z{y) to start at 
the origin and grow monotonically without limit as y — > 00. When the field is finite at 
r = (case A-2) the curve z{y) must increase monotonically from the origin and diverge at 

yo = 0'(O), showing a vertical asymptote there. In this case / (— y^) diverges at y = y^ and 
the uniqueness condition requires the lagrangian function to remain undefined for y > y^, for 
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Fig. 2. Characterization in the y-z{y) = y f (— J/^) plane of the admissible models supporting finite-energy SSS solutions. 
Observe that all curves start at y = 0, z = (y = <^ (r — ► oo) = 0, cases B-1, B-2 and B-3) and grow monotonically without 
limit. This guarantees the existence of a single cut point with the z = constant lines and thus a single- valued solution (compare 
to figure 1). 

the function z{y) to exhibit an unique growing branch. Thus in all A-2 cases the set X < —Uq 
must be excluded from the domain of definition [Vt) and then, the associated models must 



be necessarily class-2 field theories I I . The monotonicity condition for the unique branch of 



z{y) in the admissible models with finite-energy solutions takes the form 
dz 



f i-y') - 2y' 7 i-y') > 0, (2.44) 
ay 

for any y > 0. This requirement coincides with the condition (12.421) for stability which, as 
already mentioned, is fulfilled by all these admissible models. To summarize, we conclude 
that all finite- energy SSS solutions of admissible scalar field theories considered in this 
section are single-branched, stable and defined everywhere. 

To close this section let us give an expression for the energy of the SSS solutions in terms 
of the function z{y), which will be useful for the explicit calculation of the soliton energy, as 
we shall see in the examples of the next section. This expression can be obtained by taking 
into account Eqs.([23D, (12. lip and (I2.43P and reads 



y{r=0) 

As can be seen from the preceding analysis, the first term in this formula vanishes for 
admissible models with soliton solutions l£J. Once the expression of the lagrangian density 
is known, the second term gives the soliton energy directly through a quadrature. 



^^It can be shown from the analysis of the possible continuations of z(y) for y > yo that the new branches of SSS field solutions 
are pathological (non-defined everywhere, unstable, or both). 

Conversely, the conditions to be imposed on the lagrangian densities of admissible models to support finite-energy SSS 
solutions could have been directly obtained from the requirement of cancellation of the first term in 1)2. 45^ . 
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3. Some examples 



As illustrative examples of the considerations of the previous section we shall introduce 
and discuss four families of admissible models representative of the different classes analyzed 
there. 



3.1. Potential corrections to the D'Alembert Lagrangian 

The first example is given by the two-parameter family of field theories defined by la- 
grangian densities of the form 

fix) = I + AX^ (3.1) 

where A is a positive constant which gives the intensity of the self-coupling. When A = 0, 
f{X) reduces to the usual D'Alembert lagrangian density (see figure 3). The values of the 
exponent a are restricted to be irreducible ratios of two odd natural numbers (we consider 1 
as odd) a = P/Q such that 

P > a > ^) . (3.2) 




Fig. 3. Functional form of the lagrangian densities for the family of models 113.11 1 (A) for a fixed value of the parameter 
A(= 1) and several values of a(= 5/3,3, 13/3, 17/3, 25/3) (B) for a fixed value of the exponent a(= 3) and several values of 
A(= 0.01, 0.1, 1, 10, 100). The dashed line corresponds in both cases to the D'Alembert lagrangian (A = 0). 

As easily verified these restrictions allow f{X) to be defined everywhere and the associated 
energy density to be positive definite and vanishing in vacuum (admissible class-1 field 
theories). The condition a > 3/2 is imposed to ensure the convergence of the integral of 
energy at r = 0, where the SSS field strengths diverge. 

This family of models can be extended to include rational values of the exponent a = 
P/Q > 3/2 with P being an even natural number and Q an odd one. Such models are 



17 



admissible for X < if we replace A > by —A in Eg. (13. II) lU- The following considerations 
are valid for the extended family. The form of the SSS solutions is obtained from the equation 
f l2.44p . which now reads 



z{y) = I + (-l)^-Uay(^-^) = A (3.3) 



with ?/(r, A) = (f)\r,A). The function z{y) shows an unique growing branch for every value 
of the scalar charge A > and, consequently, there is an unique SSS solution of Eq. (13.31) . 
which vanishes as ~ r^^ when r — oo (case B-2 above, asymptotically coulombian) and 
diverges as ~ ^-2/{2a-i) ^]-^g^ ^ ^ q Thus the argument (X) of the lagrangian ranges 
from zero to — oo in this interval and, as expected, the stability condition (12.421) . which now 
reads 

/ (X) + 2X 7 (X) = ^ + {-lf-^\a{2a - 1)X''-^ > 0, VX < 0, (3.4) 
is fulfilled there. 

The energy of the soliton, as a function of the model parameters, can be explicitly obtained 
from the integral term of Eq. (12.451) . The final result is 

e=l^ . flf-i-.i^). (3.5) 

where B{x, y) is the Euler integral of first kind 
1 

B{x,y) = J dtf'-^l - ty-\ Re{x) > 0, Re{y) > 0. (3.6) 



In figure 4 we have plotted this energy as a function of the exponent a, with the coupling 
constant A as a parameter. We see that the energy diverges, for any value of A, as the 
exponent a approaches the value 3/2. This energy is strongly reduced in the region of values 
of a > 3/2 as the coupling constant A increases, reaching minima which vanish as A — > cxd in 
this region. When the exponent a increases the energies of the solitons become less dependent 
on A and approach asymptotically the value ej Is?!'^ = IGvr a/2/3 as a — cx). An interesting 
feature of these models is the existence of soliton solutions for any a > 3/2, no matter 
how small the coupling parameter A may be. This implies that any small correction of this 
kind to the "bare" D'Alembert lagrangian leads to the possibility of excitation of soliton 
modes. For a fixed value of the exponent, the masses of such modes increase as the intensity 
of the coupling is reduced. This behaviour is similar to the one encountered in one-space 
dimensional models supporting topological soliton solutions [32], [ID]- 
The preceding analysis can be generalized to the case of theories whose lagrangian densities 



for X < take the form 



But, for complete admissibility, this function should be matched for X > with another function satisfying the condition 
I I2.I9I I. For example, the lagrangian obtained by replacing A by A • sign{X) in 1 13.111 . Since the structure and energy of the SSS 
solutions are determined by the form of the lagrangian density for X < 0, the SSS solutions of these models are also solitons. 
^""Obviously, for X > they must be extended in agreement with the requirements of section 2. 
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Fig. 4. Energy of the solitons of the family l|3.H l as a function of the exponent a with the coupling constant A as parameter. 
The dashed line correspond to the asymptotic limit of the energy for a — > oo. The small plot shows the behaviour of the energy 
for strong self-couplings. The energy is never zero for any A < oo. 

X ^ 

MX) = - + ^(-l)^"~U„X-, (3.7) 

^ n=2 

where, for admissibility, the are constrained to form a finite sequence of positive constants 
and the exponents a„ = Pn/Qn to form an increasing sequence of rational numbers, built 
as irreducible ratios of odd natural numbers or of even and odd naturals, such that a„ > 1 
for n < N and > 3/2. In these models the energy density is positive definite for X < 
whereas the SSS field solutions diverge at the origin as 0'(r — > 0) ~ 7--2/(2a]v-i) g^j^^ vanish 
asymptotically as (p {r oo) ~ r~^. The associated energies are finite and the corresponding 
solitons are stable. 

A particular case of Eq. fl3.7p is obtained when the exponents are a finite sequence of 
consecutive naturals a„ = n{n > 1). Moreover, let us assume that we take the limit N oo 
and that the infinite sequence of A„ > converges to zero, in such a way that the series 

•y- oo 

f{X) = - + j:i-^)'"'^nX\ (3.8) 

^ n=2 

be convergent in an interval including VX < (case A-1), or in an interval including the 
range > X > — but excluding the values X < — (case A-2). We can then extend the 
family (13 .Zp to a large class of analytic functions. Such functions must satisfy the conditions 
established in section 2.2. for the integrals of energy associated to the SSS solutions of such 
models to be finite. If we assume the positivity of the coefficients A„ in Eq. (l3.8p we can 
explicitly check for these models the fulfillment of the condition of stability (I2.42p . which 
now takes the form 
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• •• 1 °° 

/ (X) + 2X / (X) = - + ^ n{2n - l)A„(-X)"-i > (3.9) 

^ n=2 

and holds in the entire domain of definition of the sohton. Thus the requirement of con- 
vergence of the series (13. 8p for any X, with the assumed restrictions A„ > (Vn) leads to 
class-1 field theories supporting SSS solitons. As an example of this let us mention the ana- 
lytic function /(X) = |s/i(X), which diverges when X — oo faster than X"' with 7 > 3/2 
(as required by Eq. (12.271) in the A-1 case) and behaves like X/2 around X = (case B-2). 

If we relax the requirement of positivity of the coefficients A„, whenever the series (13. 8p 
converges in some restricted interval X > — and the sum remains a monotonically increas- 
ing function of X there, we are lead to admissible class-2 field theories exhibiting SSS soliton 
solutions. An example of this case is the lagrangian density /(X) = ^tg[X) restricted to the 
interval — | < X < |. This model supports finite-energy stable SSS solutions belonging to 
cases A-2 and B-2. 

Let us consider the case A„ > 0,V?t,. The partial sums in (13.81) give rise to an infinite 
sequence of admissible lagrangian models of the form (13.71) with natural exponents, all of 
them supporting SSS soliton solutions belonging to cases A-1 and B-2. The explicit forms 
0Ar(r, A) of these solutions can be obtained by solving the equation (I2.43P for each lagrangian 
in the sequence. These equations take the form 

nUcl^^f^-^^^ = A (3.10) 

If the series (13.80 converges in X < 0, defining an analytic lagrangian density function 
there, the sequence of SSS soliton solutions of the partial-sum models in the expansion, 
corresponding to the same value of the scalar charge A for all N, must converge to the SSS 
soliton solution (with the same charge) associated to this lagrangian density {(j)N^oo{r, A) —>■ 
(j){r, A)). This can be directly established from Eq. (l3.10p . which defines the forms of the SSS 
solutions. Then the limit solution can be written as a functional series expansion in terms 
of the members of the sequence as 

00 

0(r,A) = 0i(r,A)+^5jv(r,A), (3.11) 

N=2 

where 5Ar(r, A) = 0Ar(r, A) — 0Ar„i(r, A). Using Eqs. (l2.1ip and (I2.12p it is easy to show that 
the sequence of energies of the equal-charge solitons associated to the partial-sum lagrangian 
densities converges towards the energy of the equal-charge soliton associated to the full series 
lagrangian (13. 8p . which can be written as the series expansion: 

00 

e(A)=ei(A)+^ Ajve(A), (3.12) 

N=2 

where AAre(A) = eAr(A) — eAr_i(A) is the difference between the energies of two consecutive 
solitons in the sequence. The first two terms of this expansion are energy- divergent. The 
first one corresponds to the divergent self-energy of the Coulomb field whereas the first 
correction A2e(A) cancels this divergence and "renormalizes" the self-energy to a finite 
value. The subsequent terms are all finite and the series converges towards the energy of 
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the limit soliton. The energy associated with the soliton of order in the sequence can be 
obtained making use of the expression ( ]2.45p and reads 

oo 

e.iA) = ^v^A^/^ / , (3.13) 



which can be numerically calculated once the coefficients are fixed. 
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Fig. 5. Energy of the sequence of SSS soliton solutions (with the same scalar charge, A) of the partial-sum lagrangian models 
obtained from the McLaurin expansion of the lagrangian density f{X) = ^sh{X), as functions of the integer parameter 
i = ' ^ being the odd exponents of the surviving terms in the expansion. As N increases these energies approach 

asymptotically the value e(A) ~ 28.5607A^/^, corresponding to the soliton energy of the exact theory. 

To illustrate this procedure let us consider the above mentioned analytic lagrangian 
f{X) = ^sh{X). The energy associated to the soliton solutions of this model, obtained 
from Eq.([235]), is e(A) ~ 28.5607A3/2_ ^jj^g partial sums of the McLaurin expansion of this 
lagrangian function are admissible models supporting a sequence of SSS soliton solutions. 
Their energies, obtained from Eq. fl3.13l) . are plotted in figure 5 as functions of i = ^y^, 
for the same value of the scalar charge A of each solution (A^ = 2i + 1 being the expo- 
nents of the surviving terms in the expansion which, in this case, are the sequence of odd 
naturals). Obviously the energy of the first-order term (A^ = l,i = 0), which corresponds 
to the Coulomb field, diverges but, as expected, the first correction already "renormalizes" 
this coulombian divergent energy and the subsequent orders reduce the (now finite) energy, 
which approaches asymptotically the energy of the soliton of the exact model as i increases. 
The convergence in this example is related to the analytic character of the sum (13. 8p but 
the "renormalization" of the divergent self-energy is due to the first correction to the pure 
D'Alembert lagrangian and would arise even if the series were not convergent. 

These results can be useful in the analysis of particle-like solutions in effective models of 
gauge-invariant interactions. Indeed, effective lagrangians arise frequently in perturbative 
schemas which lead to polynomial expressions in certain field invariants. For example, in the 
case of QED the perturbative expansion of the photon effective action, which is obtained by 
integrating out the high-energy degrees of freedom of the electron sector, defines a sequence 
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of lagrangians which take this polynomial form in the field invariants (Euler-Heisenberg 
lagrangians [32] and the higher-order corrections [33]). On the other hand, as we shall estab- 
lish in the following sections, the solution of the electrostatic spherically symmetric problem 
for a generalized gauge-invariant lagrangian model can be reduced to that of an associated 
scalar field model, whose lagrangian density is univocally defined from the gauge-invariant 
one. If the sequence of gauge-invariant effective lagrangians are of polynomial forms in the 
field invariants, the associated scalar lagrangians are also polynomials in the kinetic term, 
taking the form of partial sums of a series as (13. 8p . In this way we have established that 
the sequence of effective lagrangians describing the low-energy photon-photon interaction in 
QED support electrostatic point-like finite-energy solutions [H]. 

Let us give a tentative physical interpretation of these results. The non-linear terms in 
these effective lagrangians describe, at a classical level, a self-interaction of the gauge field 
mediated by the Dirac vacuum. The point-like solution of the "bare" Maxwell lagrangian is 
the Coulomb field, which has a divergent self-energy. The first non-linear correction term of 
the effective lagrangian (Euler-Heinsenberg) incorporates polarization effects of the vacuum 
on the classical field of the point charge, calculated to lowest order in a perturbative ex- 
pansion. These screening effects "renormalize" the charge field, which becomes finite-energy. 
The subsequent corrections in the expansion describe higher-order approximations to the 
behaviour of the screening, but the finite-energy character of the screened fields is preserved 
to all orders. Unfortunately, the validity of this effective approach is limited to energies much 
lower than the electron mass [HT] and is not accurate to describe the strong fields arising 
near the center of the particle-like solutions. Consequently, this tentative interpretation can 
not be maintained only on these grounds. For a more rigorous investigation this question 
should be considered starting from a different effective approach incorporating the vacuum 
polarization effects in presence of the strong fields of point-like charges. The analysis of this 
approach is currently in progress (see also the comments on this point in the last section). 



3.2. B-I-like models 

The second example is a two-parameter family of field theories defined by lagrangian 
densities of the form 

where is a real constant. The admissibility conditions require the values of the parameter 
a to be restricted to the range 1/2 < a < 1. Indeed, ifl<a<3/2 the energy of the SSS 
solutions diverges around r = 0. On the other hand if a > 3/2 the solution (r) is multi- 
valued and the different branches are either unstable or energy- divergent. We also discard 
the models with < a < 1/2, since the energy density in this case is not positive definite for 
X > 0. Moreover, if a is a rational number built as the irreducible ratio of an even natural 

and an odd natural numbers, the function f{X) is defined everywhere, but / (X) changes 
sign in X = —l/fi"^ and the energy becomes negative for large negative values of X. Finally, 
if a is the irreducible ratio of two odd naturals we are lead to models which exhibit multi- 
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branched SSS solutions U£J. Consequently, we must exclude the models with these values 
of the parameter a and restrict the family to the lagrangian densities which result from 
irrational values of a or rational values which are irreducible ratios of an odd and an even 
natural (always within the range 1/2 < a < 1). In these cases the lagrangian densities are 
defined only for X > — and behave as X/2 around X = 0, corresponding to class-2 
field theories and cases A-2 and B-1. There are now unique, stable and finite-energy SSS 
solutions for each model, which are defined everywhere and fall inside the cases A-2 and B-2 
(with the maximum field strength (0) = and coulombian behaviour at infinity). In the 
limits /i — or a — i> 1 Eq. fl3.14l) reduces to the D'Alembert lagrangian density. The scalar 
Born-Infeld model is a member of this family, corresponding to the frontier value a = 1/2. 

In calculating the energy of these soliton solutions as a function of the model parameters 
we evaluate the integral in 02.451) . as in the previous example. The final expression is 



3 |^|i/2 V2' 2 

We see that, as a function of fi the soliton energy behaves like e ~ and diverges as fi —>■ 

(D'Alembert limit), whereas it vanishes in the strong-coupling limit fi —>■ oo. This energy is 
not very sensitive to the exponent a in the range of admissible values and the behaviour for 
the whole family is similar to that of the scalar Born-Infeld model. 

3.3. A three-parameter family 

The third example is the three-parameter family of models defined by the lagrangian 
densities 

where a is chosen as the irreducible ratio of two positive odd naturals (in order for the 
lagrangian to be well defined on both sides of X = 0). The exponent (3 must be chosen 
as a positive irrational number or as the irreducible ratio of an odd and an even natural 
numbers. In this way the lagrangian is defined only for X > — l//i^, thus avoiding a non- 
positive definite character of the energy as well as a singularity inside the domain of definition 
These restrictions lead to a family of class-2 field theories belonging to the A-2 case. 
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being examples of the sub-cases whose lagrangians diverge at X — — l//i^ in the boundary 
of the domain of definition. The behaviours of the lagrangians around X = (r — oo for 
the SSS solutions) belong respectively to the cases B-1 (a > 1), B-2 (a = 1) or B-3 (a < 1), 
corresponding to asymptotic dampings of the soliton field strengths which are slower than 
coulombian, coulombian or faster than coulombian, respectively (see figure 6). 



^''In fact, the function z{y) in 112.441 1 has now two separated branches. The field associated to the first branch ranges in the 
interval < </> (^) < I/a* a-nd satisfies the condition 12.421 1 there, leading to a stable and finite-energy SSS solution, finite and 
defined everywhere (these branches fall inside the cases A-2 and B-2, with (p (0) = l/fj, and coulombian asymptotic behaviour). 
The remaining branch of z{y) ranges in the interval l//i < j/ < oo and exhibits a minimum at y = {^\/2a — 1)""^. Consequently 
there are two additional solutions cf) (r) defined only inside the interval < r < V2/xA(2 — 2a)('^~"'/^ (2a — l)(2o'-i)/4, 
^'^We emphasize that, as mentioned in section 2, we regard as acceptable singularities of the lagrangian density only those lying 
on the boundary of the (open and connected) domain of definition. In fact, one of the motivations in introducing this example 
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Fig. 6. Form of the lagrangian functions corresponding to models of the family 13.1611 for three sets of values of the parameter 
a ^= . All lagrangians diverge at X — > — which corresponds to the maximum value of the field strength. At the point 

X = 0, which determines the asymptotic behaviour of the solitons, the three sets of values of a give the three different 
behaviours: case B-1 {a > 0), case B-2 {a = 0) and case B-3 (a < 0). The dashed lines show the behaviour of the lagrangian 
function at large positive values of X for the admissible models corresponding to different relations among the parameters. 

We also impose the condition a > /? + 1/2, necessary to ensure the positivity of the energy 
density for any X & Q. Moreover, the convergence of the integral of energy for the SSS 
solutions as r ^ oo requires a < 3/2, as can be easily verified from the analysis of the field 
equation (12.41) and the integral of energy in this limit. Let us summarize in the following 
equations the restrictions imposed on the parameters of the models (I3.16P in order to obtain 
admissible models with soliton solutions: 

1 odd 3 odd . . , n 1 /o i ^,\ 
- < a = — — < - ; p = or irrational ; 0<p<a <1. (<J-1' ) 

2 odd 2 even 2 

where the terms "odd" and "even" are implicitly understood to apply for natural numbers. 
The D'Alembert lagrangian is a limit member of this family obtained as a — > 1 and /5 ^ 
or as a ^ 1 and /i — > 0. 
As results from the analysis of the A-2 cases, near the center the SSS solutions behave as 

0'(r ^ 0) ~ - - Xr", (3.18) 
/i 

where the exponent a is given by 

1< a = ^ < 2, (3.19) 
and A is a positive constant, which is the solution of the equation 



is to show how models with such a kind of singularities can also lead to physically reasonable results. 
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2a/xA + /3 = A/i("+^)(2A) 



(3.20) 



As in the preceding examples, the energy of the sohton solutions can be explicitly obtained 
from (I2.45p . The final expression is 



4v^7rA3/2 



a 



1/2 



^^^^ 3-2aP + 3\ 



1 3 -2a 9-2(a-/5) a - 



a 



(3.21) 



where B{x, y) is again the Euler beta-function and Ff{a, b, c, z) is the hypergeometric func- 
tion defined as 

1 } 

Fl = F{a, b, c, z) = — / t''-\l - ty-^-\l - tz)-^dt, Re(6) > 0, Re(c) > 0.(3.22) 

B[b, c-b) 

In figure 7 we have plotted the energy as a function of a with /i and (3 as parameters. For 
a given value of /i the energy is rather insensitive to the parameter (3. As a approaches 
the limit 3/2 the energy diverges for all values of /x and (3. As a function of /i the energy 
decreases as the power 1/ ji^2~°'\ 




Fig. 7. Behaviour of the energy as a function a for five values of fi and three of /3(0.01, 0.4, 0.7). Note that the lower branches 
for each value of fj. do not cover all the range 1/2 < a < 3/2 since the constraint a > /3 + 1/2 must be always fulfilled for 
admissibility. 



3.4. Short-ranged solutions 



As a fourth example let us look for a family of models whose SSS solutions are exponen- 
tially damped for large r [IfJ. In obtaining these models we shall proceed backwards, looking 



As already mentioned, these kinds of theories allow to describe short-range interactions through the exchange of self-coupled 
scalar fields. Their generalizations to the case of gauge fields, performed in section 5, lead to effective non-linear lagrangians 
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for lagrangians whose associated field equations have prescribed SSS solutions. In this way 
we shall look for a family of lagrangian density functions of the form ( 12. ip whose associated 
SSS field solutions have the simple exponentially damped formP^ 

0'(r,A)=Aexp(-a-^), (3.23) 

where A, a and n are positive constants determining the different models within this family 
P^ . The constant A is the integration constant of the first-integral (12.41) of the field equation 
(whose solution is required to be (I3.23P ) and parameterizes all SSS solutions of a given 
model. It is explicitly introduced in (13.231) by implementing the scale law (12.71) . These fields 
belong to the cases A-2 (finite field strength at the center) and B-3 (asymptotic damping 
faster than coulombian), and the corresponding models are class-2 field theories. 

By eliminating r between (I3.23P and the first-integral (12. 4p we obtain the form of the first 
derivative of the lagrangian density 

• A (2a-) 

which holds in the interval —A^ < X < (where the SSS solution (13.230 is defined) and 
diverges at the boundaries. The lagrangian density in this interval is 



/(X) = 2a2/" / = -2Aa2/" / ^dz. (3.25) 

^ ' J In2/"(A/|/) J z^l^ ^ ' 

ln(A/v^) 



As easily verified, /(O) = for any set of positive values of the parameters. In the lower 
boundary of the interval we have /(X = -A^) = -2y4cr^/"r (l - |) for > 2 (with T(t) = 
z^~^e~^dz,t > being the usual Euler gamma function) and /(X —A"^) — oo for 
n < 2 (see figure 8). This expression of the lagrangian density could be continued to the 
region X > by matching (I3.25P to any function satisfying the admissibility conditions 
there, but such continuations do not affect the structure of the solitons, which is completely 
determined by the part (13.250 of the lagrangian density 
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We calculate now the energy of these soliton solutions starting by convenience from 
Eq. (l2.10p (although formula (12.450 would also work). After a partial integration we obtain 



An 



oo 



-r'f (-<p'\r)) |^+ / r^^dr 

(Jjl 



(3.26) 



which also describe short-range interactions and preserve the explicit gauge-invariance. From this point of view such models 
provide alternatives to the usual symmetry breaking mechanism in the description of weak interactions. 

More complex choices of exponentially damped SSS fields (as, for example, if> (r) = a(r) exp (— cr j^„/2 ) i where a(r) is assumed 
to be a bounded function) may be analyzed in a similar way, but in the present example the calculations can be performed in 
terms of elementary functions. With this choice the SSS field will be a soliton, but the method works also in obtaining models 
with exponentially-damped SSS solutions which are energy-divergent. 

The constant j4 is a parameter of the model and not an integration constant of the solutions. It plays the role of the maximum 
field strength and is shared by all SSS solutions of a given model, but differs for the various models in the family. 

Nevertheless, the requirements of positivity of the energy and vanishing vacuum energy of the complete lagrangians are still 
necessary for the stability of the solitons (see section 6). 
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A 




Fig. 8. Functional form of the lagrangian densities for the family 1)3.2511 . exhibiting short-ranged soliton solutions. The dashed 
line indicates possible continuations of the lagrangian density for X > 0. 

The integrated part in this equation can be shown to vanish for the prescribed solutions. In 
calculating the integral in the second term we use the first-integral field equation (12.41) and 
the expression of the field (13.231) . which leads to 

£ = -2f (X)0'(r)0"(r) = -^<P"{r) = 2A^'-y'Anar-''exp [-ctj^) , (3.27) 

and the final expression for the energy of the solitons reads 

SttA fl + n\ , , 



3ai/" 



n 



which is proportional to the maximum amplitude of the field strength and decreases as the 
range of the field is reduced. In figure 9 we have plotted the energy of the unit maximum- 
field-strength as a function of n for several values of the constant a. As we see the energy 
diverges as n — > and becomes less dependent on the constant a for large values of n, 

e — Stt 
AA372 3 



approaching asymptotically a fixed value " — ^'^ 



4. The multicomponent scalar field 



We shall extend the results of section 2 to the case of a set of scalar fields 0j(x)(z = 
1...X). In many cases the covariant lagrangians including N scalar fields and their first- 
order derivatives are constrained by conditions which allow to implement some internal 
symmetries. Such conditions manifest themselves in the structure of the manifold where 
the field takes its values. A well-known example is the non-linear sigma model where this 
manifold is a Riemann space implementing chiral symmetry [n],[Sn]- Here we shall restrict 
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Fig. 9. Energy for short-ranged solitons 113.231 1 of unit maximum-field strength as a function of the exponent n, for several 
values of the parameter a (note that we have plotted ^^^3/2 instead of j^y2 ''^ previous examples). 

ourselves to the case where the field manifold is the A^-dimensional Euclidean space and the 
SO{N) invariant lagrangian density depends only on derivative terms 




5/^00=/ 2: W^^n, (4.1) 



where, as in section 2, f{X) is a given continuous, derivable (C^ for X < and monotoni- 
cally increasing > 0, VX 7^ 0; ^ > 0,X = 0) function defined in a open and connected 
domain (f2 C 3fJ) which includes the vacuum ((X = 0) G Besides the fact that these 
models are the natural generalizations of the scalar field theories studied so far, there is an- 
other motivation for their analysis. Indeed, as we shall see in the following sections, a class 
of soliton solutions arising in generalized gauge field theories of some compact semi-simple 
Lie groups of dimension N reduce to multiscalar {N components) solitons of some of the 
models (H?T]) . 

The field equations associated to the lagrangians (14.11) take the form of N local conservation 
laws 

d,J^ = 0, (4.2) 
where the conserved currents Jf are 

Jt =f {X)dy,, (4.3) 
with X = J2f=i da4'i ■ <9"0j. The canonical energy-momentum tensor is 

T^. = 2f (X) J2 d^(j)A(l)^ - /(X)r/^„ (4.4) 

i=l 

and the corresponding energy density 
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p(^) = 2/WE(^) - fix), (4.5) 

is positive definite under the same conditions established for the lagrangian function f{X) 
of the one-component case. 

For the SSS solutions 0i(r), equations (14.21) have N first-integral field equations of the 
form 

rV;/ (|-E0?j =A., (4.6) 

where (f)[ = d(j)i/dr and Aj are the integration constants. Now the signs of every component 
of the scalar field and of the corresponding integration constant are the same, but may differ 
for different components. In order to solve the system (14. 6 p let us introduce the functions 
Xi{r) = — 0^^(r), in such a way that X{r) = Y,f=iXi{r). By squaring and adding Eqs. (l4.6p 
we obtain 

,2 N 

r'Xf {X) = -Y^Al (4.7) 

i=l 

or, equivalently. 



r'^rXf{X) = A, (4.8) 
where 



A 



N 



EA?- (4.9) 



Equation (14. 8 p has the same form as the first-integral of the one-component scalar case (12.41) . 
Consequently, if the function f{X) is the same in both cases, we can associate to any SSS 
solution of the one-component case, of the form (p (r. A), a set of sequences of functions 
which are SSS solutions of the multicomponent scalar equations. Such functions take the 
form 

0;(r,A,) = ^0'(r,A), (4.10) 

and, owing to Eq. (l4.9p . there is a one-to-one correspondence between such sequences and the 
points of the sphere of radius A in the A^-dimensional Euclidean space (Jft^). Obviously, this 
is a straightforward consequence of the invariance of the lagrangian (14. ip under rotations 
in the internal space. The constants A, can now be identified as the "source point-charges" 
associated to the different components of the SSS field, namely 

K = ^ I d'xV ■(} {X)Vct>^. (4.11) 



An 

The field potentials obtained by integration of (I4.10p read 

0,(r,A,) = ^0(r,A) + A„ (4.12) 
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where Aj are integration constants and A is the mean-square scalar charge (14. 9p . 
If we consider now the energy associated to these SSS solutions we find from ( 14. Sp 

e= -An J I ^Xi(r) Mr = -An J r"^ f {X{r)) dr. (4.13) 
^*=i ' 

This is the energy of the one-component SSS solution corresponding to the integration 
constant (14.90 . Thus the set of SSS solutions of the multicomponent scalar field associated 
to the points of the sphere of radius A in 3?^ is degenerate in energy. Moreover, the search 
for conditions to be imposed on f{X) for the existence of finite-energy SSS solutions (as 
well as the admissibility constraints) in the multicomponent case reduces to the analysis of 
section 2 for the one-component case. 

Concerning the conditions for stability of the solutions (I4.12p . the analysis of the one- 
component case can be straightforwardly generalized to the present situation (see subsec- 
tion 6.3. for details). The final conclusion is that the multicomponent soliton solutions of 
admissible models are also linearly stable against charge-preserving perturbations. 



5. Gauge fields 



We shall now extend the analysis developed in the previous sections to generalized gauge 
field theories of compact semi-simple Lie groups. To start with, we study the simpler case 
of generalized ([/(l)-invariant) electromagnetic fields before proceeding further to the case 
of non-abelian gauge fields. 



5.1. Abelian case 



We define lagrangian densities for generalized electromagnetic fields defined as arbitrary 
functions of the two quadratic field invariants, built from the Maxwell tensor and its dual. 
Following the conventions of Ref. |12], these tensors are defined as 



F^, = d^A,-d,A^ (5.1) 



where _ —^^-^^^ = 1. The electric and magnetic fields are defined as E"^ = — and 
= —^e'^^^Fjk while the quadratic invariants X and Y are 

X = -^F^,F^'' = E^-H^ (5.2) 

Y = ~F^,F*f''' = 2E-H. 

We define the general form of the lagrangian density as 

L = ^{X,Y), (5.3) 

where <y9 is a given continuous and derivable function on its domain of definition (Q) of the 
X — Y plane (3?^). As in the scalar case we assume fl to be open and connected and including 
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the vacuum ((X = 0,Y = 0) E Q) . As a. minimal extension of the assumptions of section 2 
concerning the regularity properties of the scalar lagrangian functions, which is necessary for 
future purposes (see Eg. (15.201) below), we shall assume ip{X, Y) to be of class on the line 
(X > 0, y = 0) n ^ and dip/dX to be strictly positive there. By generalizing the definitions 
of section 2 we shall call "class-1 field theories" the models defined and regular everywhere 
{^l = 3?^) and "class-2 field theories" those with C 3?^. We also require (p{X, Y) to be 
symmetric in the second argument, in order to implement parity invariance 

^{X,Y)=^{X,-Y). (5.4) 

This implies that the odd partial derivatives of (p{X, Y) with respect to Y must vanish on 
Y = 0. The Maxwell lagrangian density corresponds to (p{X, Y) = while the Born-Infeld 
electrodynamics is given by the lagrangian density 
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Lb-i = fB^iiX, Y) = i— , (5.5) 

where ^ is the maximum field strength, attained at the center of the solution. The lagrangian 
(15. 5p reduces to the Maxwell one in the limit /i — > 0. 

The symmetric (gauge-invariant) energy-momentum tensor obtained from the lagrangian 
density (15. 3p takes the form 

= 2 f||^M«^" + I^^M^^r) - = 2||^M«^" + - v] V,u, (5.6) 



and the associated energy density is 

/ = IS. = 2||b- + 2p . H - ,(.Y. Y) = 2.y|| - ^(A-. Y) + + 2||//^(5.7) 

We assume the symmetric energy-momentum tensor (15. 6p to give the correct space-time 
energy density distribution. Let us analyze the conditions for the positivity of the energy 
density of any field configuration. The inspection of Eq. (l5.7p . together with the requirement 
of vanishing of the vacuum energy, lead to the set of necessary conditions 

<^(0,0)=0 ; (p{X,0) <0 y{X <0,Y = 0) en ; 1^ > V(X,r) G fi, (5.8) 

to be satisfied by the lagrangian densities. However, it is possible to obtain a minimal set 
of necessary and sufficient conditions of admissibility for a satisfactory energetic behaviour. 
Solving Eqs. (l5.2p for the fields we obtain 



Y2 




+ x] > UVX^+Y^ + X) (5.9) 



cos2(^9) - 2 



1 



X\ > - (VX2 + Y^-X 

where is the angle between E and H. From these expressions the energy density can be 
written as 
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Consequently, the requirement of the positive definite character of the energy leads to the 
minimal necessary and sufficient condition 



> (Vx^TT^ + x) ll + y 1^ - ^(x, y) > 



dX 



dY 



(5.11) 



to be satisfied in the entire domain of definition {Vt). Generalizing the criteria of section 2, 
we only regard as "admissible" those models whose lagrangian densities satisfy the condi- 
tion (15. lip , aside from the vanishing of the vacuum energy and the regularity and parity- 



invariance conditions stated above I I . Consequently, the admissible lagrangians must be 



solutions of the first-order linear inhomogeneous partial differential equation 

dX dY 



(VX^TT^ + X)^ + Y^- ^(X, Y) = v^(X, Y), 



(5.12) 



where \1/(X, Y) is any function being positive definite in Q and vanishing in vacuum (^'(O, 0) = 
0). Such solutions must also satisfy the regularity and parity- invariance requirements, as sup- 
plementary conditions. 
Let us now consider the trace of the symmetric energy-momentum tensor 



dip 
dX 



{^E^^H^)+2^E-H-ip{X,Y) 



.(5.13) 



From the last expression we see that the sub-class of models with traceless symmetric energy- 
momentum tensors is given by the lagrangian densities 'p{X, Y) which are solutions of the 
first-order linear homogeneous partial differential equation 



dX dY ' > 



0. 



(5.14) 



The general solution of this equation is the family of all conic surfaces in the (X, F, y9)-space 
having the origin as a vertex. Clearly the set of planes of the form ip = aX + bY (a and b 
being constants) are particular solutions of this equation (which violate parity invariance if 
6 7^ 0). The simple case 6 = with a = ^ corresponds to the Maxwell theory. 
The field equations obtained from the lagrangian (15. 3p are 



"""[dX ^ dY , 



0, 



to which we add the Bianchi identities 
0. 

In terms of the fields these equations read 



d^F*^"' 



(5.15) 



(5.16) 



^Obviously, excepting the vanishing of the vacuum energy, the remaining conditions in 1 15.81 1 are consequences of (15.111 1. 
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and 

^ dH 

V X E-- 



dt 

V-H = 0. (5.18) 

For electrostatic fields we have Y = and in the ESS case these equations can be written 
in terms of the electrostatic potential A^{r) {A = 0,E{r) = — VA'^(r)) in such a way that 
the first group of equations in fl5.17p leads to the first-integral 
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where now X = and q is an integration constant. Using the identification <^(r, A) = 

A°(r, g = A) this equation coincides with the first-integral (12.41) for the SSS solutions of a 
scalar field model with a lagrangian density defined by 



Lscaiar = f{d^(l> " 9^0) = f (X) = -v{-X, Y = 0), (5.20) 
which leads to 

fiX) = ^{-X,Y = 0). (5.21) 

Conversely we can associate to each scalar model defined by a lagrangian density f{X), 
a family of electromagnetic field models defined by lagrangian densities if{X, Y) satisfying 
Eq. (15.201) as well as the admissibility (15. lip and stabihty constraints (see Eq. (15.301) below). 
The ESS field solutions of all electromagnetic generalizations {\E{r, q)\) have the same form, 
as functions of r, as the SSS field solutions (0 (?", A)) of the original scalar model, q and A 
being the integration constants, which should be identified as the electric and scalar point-like 
charges associated to the solution, respectively. Indeed, in the generalized electromagnetic 
case the definition of the electric charge associated to a given field is 




d'rV- r^E+^H I (5,22) 



1 

Air 

which now is conserved as a consequence of the field equations. Substituting in this equation 
the ESS field coming from the solution of (I5.19P we obtain q as the value of its total electric 
charge 



4:71 



dX 



g53(r). (5.23) 



^^The remaining three equations in 115.1711 are identically satisfied for arbitrary electrostatic fields, owing to Eg. 115.41 1. On the 
other hand, the set of equations 1 15.181 1 are trivially satisfied by the ESS solutions. 
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As pointed out in Ref. [3] for the Born-Infeld model this charge can be interpreted as a 
source point-hke charge at the center of the ESS solution or, alternatively, as a continuous 
charge- density distribution associated with the field and given by 

1 dip 



.(X = o,y = 0)V-E. 



(5.24) 



AndX 

This interpretation, as already discussed for scalar models, requires the function r'^E{r) to 
vanish at the origin (this condition is always fulfilled for the finite-energy ESS solutions) 
and the field E{r) to be asymptotically coulombian (B-2 case models). 

In calculating the total energy of these electrostatic central fields from the energy density 
(15. 7p we are lead to 



eM = Svr / r^ll \x = E\r, q),Y = o] E\r, q)dr 



4:71 / rV X = E'^{r, q),Y = dr, 



(5.25) 



(the index e stands for electric field). The energy associated with the corresponding SSS 
scalar field solutions, obtained from Eqs. (l2.10p and (I5.20p reads 



e,(A) =47r J r'if [E\r,A),0 





dr, 



(5.26) 



(the index s stands for scalar field). If the total energy fl5.26p associated to the scalar field 
is finite, so is ee(g). Indeed, using Eq. (l5.19p this energy becomes 



ee{q) = Siiq \A%oo, q) - A°(0, q)] - e,(g) 



which, owing to (I2.12p . is related with the scalar soliton energy through 
ee{q) = 2e,(g), 



(5.27) 



(5.28) 



and must be also finite. Thus the energy of a ESS solution is twofold the energy of the 
corresponding SSS soliton when the integration constants take the same value (A = q). 
Equivalently, Eq. fl2.10p . which gives the scaling of energies, leads to the following relation 
between the integration constants (the electric (g) and the scalar (A) charges) of an electro- 
static soliton and the associated scalar soliton of equal-energy 

2)2/3A. 



(5.29) 



We then conclude that the ESS solutions for the families of electromagnetic models which 
generalize (through Eq. fl5.20p ) the different classes of scalar models with soliton solutions, 
have the same functional forms as the corresponding SSS scalar solutions and are also of 
finite-energy. Moreover, the classification of the admissible models with soliton solutions in 
the scalar case according to the central and asymptotic behaviours of the fields, immediately 
induces, through Eq. fl5.2Up . a similar classification of the finite-energy ESS solutions in the 
electromagnetic case. 
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Let us point out an immediate consequence of this analysis (which is a corollary of the 
non-existence theorems established in Refs. [20j,[43]): there are not ESS soliton solutions for 
admissible generalized electromagnetic field theories with traceless energy-momentum tensor. 
Indeed, as mentioned above the lagrangian densities of such theories (see Eq.f l5.13p ) are given 
by conic surfaces in the (X, Y, ip) space and the associated scalar field lagrangian densities 
/(X) = — (/?(— X, Y = 0) are straight lines in the (X, /) plane for X < 0. Consequently the 
associated SSS solutions are coulombian in form as well as energy-divergent, and so are the 
ESS solutions of these generalized electromagnetic models. 

Although finite-energy SSS solutions of admissible scalar models are always linearly stable 
against charge-preserving perturbations, this is not so for the finite-energy ESS solutions of 
admissible generalized electromagnetic theories. Indeed, the analysis of the linear stability of 
the electrostatic solitons leads to a generalization of the criteria obtained in the scalar case 
(see section 6 for details). As a result of this analysis, the electrostatic finite-energy central 
field solutions of admissible generalized electromagnetic field models, with their lagrangian 
densities (p{X, Y) satisfying the supplementary condition 

i--|:i>0. (5.30, 

in the entire domain of existence of the ESS solutions in the plane (X, Y = 0), can be shown 
to be local minima of the energy functional against small charge-preserving perturbations. 
Consequently, admissibility, finite-energy condition of the ESS solutions aside from Eq. fl5.30p 
are necessary and sufficient conditions for linear static stability. Moreover, the linear analysis 
of dynamics of the small perturbations of the soliton solutions performed for the scalar 
models can be generalized to the associated families of electromagnetic models which satisfy 
fl5.30p (see section 6). 

Finally, the conditions for univoque and everywhere defined ESS solutions are straightfor- 
wardly deduced from those of the scalar case, aside from Eq. fl5.20p . 

5.2. Non-abelian case 

The results obtained for generalized electromagnetic fields can be extended to non-abelian 
generalized gauge field theories of compact semi-simple Lie groups of dimension N. As 
usually, in this case the tensor field strength components in the algebra and their duals are 
defined from the gauge fields Aa^ and the structure constants Cabc as 

be 

Kixu = 2^tMuaf3F°'^, (5.31) 

whose components define the fields Ea, Ha in the usual form. In order to introduce the 
lagrangian densities governing the generalized dynamics of these fields we must define per- 
tinent field invariants. However there is now an ambiguity in the calculation of the traces 
over the group indices, leading to different possibilities in the definition of these invariants. 
Although at this regard different prescriptions have been introduced, mainly in the context 
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of string theory, where B-I-hke actions arise as a low-energy effective field limit [16], [H], here 
we shall restrict our analysis to the case of actions built from the two simplest first-order 
field invariants, defined from the ordinary prescription for the calculation of the traces as 

^ = -^E(Wr)=E(^a-^a) 
^ a a 

Y = -l Y.iFa,uFr) = 2 E (^a ■ Ha) , (5.32) 

where 1 < a < A^. The generalized lagrangian density is now assumed to be a given function 
93 (X, Y) which (again for parity invariance) must be symmetric in the second argument 
((/j(X, Y) = ip{X, —Y)) and satisfy the same admissibility constraints of definition, continuity 
and derivability (as well as the distinction between class-1 and class-2 field theories) as in 
the electromagnetic case. 
The associated field equations read now 



E^^ 



dX ' dY ^ 



0, (5.33) 



where Dacfi = Sacdfj, — gJ2b Cabc^b^ is the gauge-covariant derivative. 
The symmetric energy-momentum tensor is 

T^u = 2 E Fa,a {^^F:^ + I^F;-^ - y,r/,„ (5.34) 

and the energy density takes the form 

/ = = 2|| E El + 2^j:Ea-Ha- viX, Y). (5.35) 

The admissibility conditions to be imposed on the lagrangian density, for the energy func- 
tional to be positive definite and vanishing in vacuum, take the same form as in the elec- 
tromagnetic case (see Eqs. (15.81) and (15. lip ). Also the trace of the energy-momentum tensor 
has the same expression (see the last member of Eq. (l5.13l) ). and vanishes under the same 
conditions (15.141) . The subclass of models with non- vanishing trace energy- momentum tensor 
breaks the scale invariance and thus circumvent the non-existence theorems [20j , [43j allowing 
for soliton-like solutions. 
Let us consider the ESS solutions of these models. We consider fields of the form 



4(r) = -V(A°(r)) =-<°(r)- ; #, = 0, (5.36) 
where the functions A^^{r) are the time-like components of the gauge potential in the Lorentz 

dr 



gauge {Aa = 0) and A'^ = When replaced in the field equations (15.331) we are lead to 




-V ■ ( ^VA'Jr)] = -V ■ ( ^Ay(r)- ] = 





dX 

.^E^a.X(r)E, = -^^a,X(r)VA°(r) = (5.37) 
be C'A 
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DC 



where X = J2a E^- The first group of equations has a set of first-integrals of the form 



r'^K^r) = Q., (5.38) 

where the Qa are integration constants which will be identified below as "source color 
charges". With the identification (pair) = A^{r) and = Qa, these equations coincide 
with the field equations (14.61) for a multicomponent SSS scalar field theory whose lagrangian 
density is given by 

^ = / (E 9,<P- ■ = f{X) = -vi-X, Y = 0). (5.39) 

Thus the solutions of equations (15.381) are obtained from Eq. (l4.10p as 

\Ea{r,Q,)\ = A'^'{r,Q,) = ^<P'{r,Q) ; i?, = 0, (5.40) 



where Q = yJ2aQl is the mean-square color charge and (j)'{r,Q) is the SSS solution of the 
associated one-component scalar model defined by a lagrangian density of the form (15.391) . 
These equations can be integrated once, leading to 

A'a{r,Qt) = ^Hr,Q)+Xa, (5.41) 

where Xa are integration constants. These functions must also satisfy the second set of 
equations (I5.38p . Owing to the antisymmetry of the structure constants these equations 
lead to the supplementary restriction 

Xa = ^X, (5.42) 

and the final solution of (15.381) is 

Alir,Q,) = ^i<Pir,Q) + x), (5.43) 

where x is an arbitrary constant, being now the same for all components of the potential. 
In terms of the fields the final solution is still given by Eq. (15.401) . The integration constants 
Qa in this solution must be interpreted as source point-like color charges associated with 
the different components of the gauge field. Indeed, color charges are defined in general as 

Q. = ^fd,fV-(%K + ^HX (5.44) 



4:71 J " \dX BY 

which, owing to the field equations (15.331) . include now external source charges and charges 
carried by the field itself. The latter ones come from the integration of the term 



9 E ^abcAbfjL 
be 



dX ' dY ' 



(5.45) 
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which vanishes for the ESS solutions. On the other hand the former are associated to Dirac 
distributions of weight 47rQa, as can be easily seen from the substitution of Eqs.( 15.38|) in the 
first set of Egs. dHIaSjl . 

The calculation of the energy associated to these solutions proceeds in the same way as in 
the electromagnetic case. The integration of Eg. (15.351) gives 

b 

(5.46) 

where the index gf stands for gauge field. From this expression it is now straightforward to 
show that, as in the corresponding multiscalar case, this energy is finite if the energy of the 
associated scalar solitons is finite, and depends on the charges only through the constant Q, 
having the same kind of degeneration on spheres of radius Q in the A^-dimensional color- 
charge space. The relation between the finite energies of the solitons with equal mean-square 
charges in the gauge models and in the associated multiscalar models is given by the same 
equation (15.281) relating the finite energies in the cases of the abelian models and their 
associated scalar models. 

Concerning the stability of the gauge field solitons we shall show in the next section that 
the finite-energy ESS solutions of admissible generalized non-abelian models are linearly 
stable if (and only if) the lagrangian density functions satisfy the same criterion (I5.30p 
obtained in the abelian case. 

Let us summarize the main conclusions of this section: The set of generalized gauge field 
theories of compact semi-simple Lie groups, whose lagrangian densities are functions (p{X, Y) 
of the field invariants (15.321) . satisfying the admissibility conditions and the stability criterion 
(I5.30p . and supporting finite-energy ESS non-topological soliton solutions, can be split in 
equivalence classes. Two models belong to the same class if their respective lagrangian 
densities satisfy the condition ipi{X,0) = ip2{X,0). The forms of the ESS soliton solutions 
and their energies coincide for all models belonging to the same class. There is a one-to- 
one correspondence, given by Eq. (15.201) . between the set of these classes and the set of 
admissible scalar field models defined by Eq. (l2.1l) and supporting finite-energy SSS non- 
topological solitons. The forms and energies of the gauge solitons are obtained from those 
of the corresponding scalar solitons through Eqs. (l5.40p . (I5.4ip and (I5.28p . The analysis and 
classification of scalar solitons performed in section 2 can be immediately generalized to 
gauge solitons through this correspondence. Furthermore, the explicit examples supporting 
scalar solitons, introduced in section 3, can also be extended to the gauge field case simply 
by including the Y invariant in such a way that the admissibility (15. lip and stability (I5.30p 
constraints be fulfilled by the extended models. 

6. Stability 

Let us now analyze the stability of the finite-energy solutions of the different models in- 
troduced so far. We make a distinction between "strong" stability, defined as the ability of 



X = Y^El{r,K),Y = Q 



4:71 / r'^ifi 



X = Y^El{r,K),Y = Q 



d\ 
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a soliton to maintain its identity under any perturbation or in closed many-soliton config- 
urations [4], [45], and "weak" stability, identified with the usual linear stability under small 
perturbations. Rigorous analysis of stability in the strong sense has been performed for a few 
field theories in one-space dimension which exhibit conserved discrete topological charges 
associated with the soliton solutions. In three-space dimensions similar topological conserva- 
tion laws are responsible for the stability of the 't Hooft-Polyakov monopole solution [5] or 
the chiral soliton solution of Deser et al. [9]. But satisfactory general methods for the analy- 
sis of interactions between non-topological solitons and strong external fields in three space 
dimensions are still lacking and only numerical analysis of the evolution of the solutions 
can give some insight on this issue for most models. In our context, a tentative approach to 
this question has been developed by Chernitsky for the Born-Infeld model [46]. It is based 
on the use of the discontinuity of the field strength at the center of static B-I solitons as 
a marker of the presence and location of the dynamic soliton evolving in interaction with 
strong external fields, or in many-soliton configurations. Since all SSS soliton solutions of 
the models considered here exhibit similar central field singularities, this procedure might 
be extended to these cases, but such an extension lies beyond the scope of the present work. 

In the case of interactions between solitons and weak external fields (or for widely separated 
soliton configurations) linear stability ensures the identity preservation of the solitons and 
becomes a basic condition for the consistency of the low- energy analysis. The results of this 
analysis may be interpreted in terms of particle-field (or particle-particle) force laws and 
describe the radiative behaviour in these processes |46j . 



6.1. Static stability of scalar solitons. 



We shall begin with the study of the static linear stability for the soliton solutions of the 
scalar models of section 2, by analyzing the behaviour of their energy, which must be a min- 



imum against appropriate small perturbations I I. We shall show that, for these models, the 



conditions of admissibility guarantee this kind of stability for all finite-energy SSS solutions. 
We start with the SSS potential 0(r) and introduce a set of small static perturbations 5(f){r), 
finite and regular (as well as their first order spatial derivatives) everywhere. We also require 
the perturbation to leave unchanged the scalar charge associated to the solution 1^ . To first 
order in the perturbations the modification of this scalar charge, obtained from Eq. (12.51) . 
reads 



An 



d^rV 



f (Xo)(V(50) - 2 7 (Xo) (V0 ■ V0 



0, 



(6.1) 



where now Xq = — (V(0))^ = — 0^(r). The condition AA - 
behaviour of the admissible perturbations at r = and as r 
satisfy 



imposes restrictions on the 
> oo. In particular, 6(j){r) must 



^*As emphasized in Ref. |47j this criterion is a sufficient condition for stability, but is not a necessary one. Here we shall not 
consider the problem of stability of solutions which do not correspond to minima of the energy, a complicated task which 
deserves a study in itself. 

This condition is essential for the energy of the soliton to be a minimum. Indeed, without such a condition the perturbation 
does not necessarily lead to an increase of the energy of the soliton (as can be easily seen by differentiating Eg. 112.101 1 with 
respect to small variations of A) and the perturbed soliton might evolve towards less energetic states. 
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lim 

r— »oo 



0. 



oo - 



(6.2) 



In this manner the perturbed fields remain inside the space of functions defined by the 
prescribed boundary conditions (on 6*00 in this case) which determine uniquely the solution 
associated to a given value of the charge. At the center of the soliton 50(r) must be regular 
(see Eg. (16.271) and the analysis of the dynamic stability below). 
The first-order perturbation of the energy, calculated by expanding (I2.10p becomes 



Aie = 2 
= 2 



rf3r/(Xo)V(0)-V((50) = 

dsfV ■ (f (Xo)50V(0)') - 2 / dsfScl>V ■ If (Xo)V(0) 



(6.3) 



where a partial integration has been performed. Owing to Eq. (12.41) and the assumed asymp- 
totic behaviour of the perturbation 50(r), the two integrals in the last equation converge 
and cancel each other so that the first variation of the energy vanishes. This is the necessary 
condition for the energy of the soliton to be an extremum. The second variation reads 

, 2 



Aae = Jd^ff {Xo){V6<pf - 2 J d^fj (Xq) (V0 ■ 



(6.4) 



where, owing to the boundary behaviours of the perturbation and the SSS field itself, both 

integrals are also convergent. From the arbitrariness of 50, the positivity of / {X) and the 
minimum condition of the energy > 0, we see that static stability is reached if the 
requirement 



/ (X) < 0, 

is fulfilled in all the range of values of X = Xq 
if we rewrite Eq.f l6.4p as 



dsf 



f (Xo) + 2X0 / (Xo) 



+ 



/ 4r/(Xo) 



^ ["oe 



dr 
1 



+ 



(6.5) 

(r) covered by the solution. However, 



(6.6) 



COS^(^^) \ (9(y9 



we are lead to the less restrictive static stability criterion 
f (Xo) + 2X0 7 (Xo) > 0, 
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(6.7) 



which is a necessary and sufficient condition for linear stability, as opposed to Eq.f l6.5p which 
is only a sufficient one. This criterion is always fulfilled for admissible models with finite- 
energy SSS solutions. Indeed, by deriving the first-integral equation (12. 4p with respect to r 
we obtain 

2A 



f{Xo 



(r) f (Xo) 



(6.8) 



'The only model for which 1 16.71 1 vanishes everywhere is singular and corresponds to the lagrangian f{X) = \\/~X 
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which, owing to the monotonicity of (j)' (r), is positive in all the range of values of Xq covered 
by the solution. We conclude that the finite-energy SSS solutions of admissible scalar 
models are statically stable. 

Against perturbations which modify the charge the solitons are unstable, but these insta- 
bilities are blocked if charge conservation is implicit in the model (as in the case of generalized 
gauge field theories considered below) or if it is a consequence of the nature of the external 
sources. 



6.2. Dynamic stability of scalar solitons. 



Let us now consider the dynamic stability of the SSS solutions. The initial perturbation 
defines the following Cauchy conditions 

(9$ 



= 0) = 0(r) + 



dt 



,(f,t = 0) =0, 



(6.9) 



for a dynamical problem determining the temporal evolution of the perturbed field ^(r, t), 
which is governed by the hyperbolic field equations (12.21) . At the first order, the evolution 
of the perturbation 



50(f, t) = $(r, t) - (j){r), 
is given by the linearized scalar field equation 
d 



dt 



dt 



V- 



/ (Xo) V(50) - 2 / (Xo) (V0 ■ V0 



0. 



This is the Euler-Lagrange equation associated with the lagrangian density 



/ {Xo)d^6<j) ■ d^5(t) - 2 7 (Xo) (V0 ■ V50 



(6.10) 



(6.11) 



(6.12) 



which is defined everywhere. Equation (16. lip has the form of a local conservation law for a 
charge density rj =f (Xq)^ which, in integral form, becomes 



^~dr 



f (Xo) V(50) - 2 / (Xo) ( V0 ■ V(50) V0 



0. (6.13) 



The r.h.s. of this equation is proportional to the first-order perturbation of the scalar charge 
of the soliton (see Eq. fl6.ip ) which, consequently, remains conserved as time evolves. More- 
over, for solutions satisfying the initial conditions (16. 9p the quantity 



rf3r/(Xo)(50(f,t), 



(6.14) 



remains constant in time. 

Equation (16. lip , together with the conditions (16. 9p and (16. 2p . outline a spectral problem 
to which we can apply standard methods. We look for solutions separating time and spatial 
variables under the form 



(5</.(f,t,r)=T(t,r)^(r,r) 



(6.15) 
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where F is the separation constant and the eigenfunction ip{f, T) is assumed to satisfy the 
boundary condition ( 16. 2p . Replacing this expression in (16.1ip and using the second of the 
initial conditions (16. 9p we are lead to 



T(t,r) 



cos 



iVft), 



and 



r / (Xo)^ 



/ (Xo)V(^) - 2 / (Xo) 



k'2 



r r 



(6.16) 



(6.17) 



The sign of the eigenvalue F is crucial for stability. Multiplying this equation by V^(r, F) and 
integrating over all space we are lead (after an integration by parts of the right-hand-side) 

to 



/(Xo)(W)^-2 7(Xo)0'2 



(6.18) 



Owing to Eq. (16. 8p . together with the admissibility and boundary conditions, both the r.h.s. 
of this equation and the coefficient of F are finite and positive. Then so is for F, and the 
evolution is oscillatory and bounded in time P^ . Moreover, if we consider two different 
eigenvalues (Fj,i = (1,2)) and their associated eigenf unctions {ipi = ip{r,Ti),i = (1,2)) 
equation (16.171) leads to 



(F2 - Fi) J d^ff iXo)^i^2 = J 4r (^2V ■ El - ^iV • 
where 

=/ (Xo) V^, - 2 7 (Xo) (V0 ■ V^,) V0. 



2 ) 



(6.19) 



(6.20) 

After a partial integration and making use of the boundary conditions we see that the 
right-hand-side of (16.191) vanishes and thus we are lead to the orthogonality relation 



dsf f {Xo)ipii'2 



0. 



(6.21) 



These results outline a Sturm-Liouville problem for each admissible scalar model of the form 
(12. ip supporting finite-energy SSS solutions and lead to the following conclusions [H]: i) The 
analysis of the dynamics of the small oscillations around these solutions leads in all cases 
to discrete spectra of eigenvalues Fj, ii) The associated eigenfunctions are orthogonal and 
finite-norm with respect to the scalar product 



J- i 



(6.22) 



defined with the kernel / (Xq) > 0. Such functions generate a complete Hilbert space in 
which any perturbation can be expanded. 



^''Note that the r.h.s. of 1 16. 181 1 coincides with the second variation of the energy associated to the eigenfunction i/)(r', F) and 
has the same sign as F. This establishes a strict correspondence between static and dynamic stabilities of the SSS solutions. 
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On the other hand, we can now separate the spatial eigenf unctions in radial and angular 
coordinates as 



7/>(r,^,yp,r,0 = i?(r,r,Z)FK^,<^), 
where the angular components are the usual spherical harmonics satisfying 



sin '(9 



d_ 



The radial components obey to the equation 



d 



1 dR\ l{l + 1) - Tr^ 



dr \ rS" dr 



2r2 



R = 0, 



(6.23) 



(6.24) 



(6.25) 



which is obtained from fl6.17p . (16.231) and using the first-integral (12. 4p . These equations have 
the standard Sturm-Liouville form US] 



Ly + \ji{x)y = 



Ly = — 

dx 



— q{x)y with k{x) > 0, fi{x) > 0. 



(6.26) 



The asymptotic behaviour of R{r) is obtained from the asymptotic form of the admissible 
solitons (0 (r oo) ~ l/rP;p> 1) through 



d^R p dR 
h - — 

dr'^ r dr 



^ 2r2 



0, 



(6.27) 



which is a Lommel equation and can be solved in terms of Bessel functions [49j . For large r 
we can assume for the solution the asymptotic form 



R(r 



oo ~ 



(6.28) 



where Q{r) is a bounded function and, owing to the boundary condition (16. 2p . q is restricted 
to be g > p > 1. By neglecting the higher-order terms in 1/r in the resulting equation for p 
we are lead to 



0. 



(6.29) 



Owing to the positivity of V the solution of this equation is oscillatory and the asymptotic 
behaviour of the eigenf unctions is given by 



28 



COS 



Rir 



oo ~ 



(6.30) 



where x is a constant phase. This asymptotic form of the eigenfunctions makes the integral 
of the second variation of the energy (16.61) to converge in the r ^ oo limit. 



^®The value of g as a function of p and I can be explicitly obtained from recursion formulae methods )50| applied to EQ.I I6.27t . 
which also allow for the approximate determination of the eigenvalue spectrum of F. 
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To determine the behaviour of R{r) around the center of the sohton we must consider 
separately cases A-1 and A-2. Let us assume in both form 

i?(r 0) ~ a - /5r^. (6.31) 

In case A-1 (0'(r ^ 0) ~ 1/r''', < p < 1) equation fl6.25p becomes 

2^g^±l^r« + (m + 1) - Tr^) (a - (3r^) ^ 0. (6.32) 
p \ ' 



For the first surface spherical harmonic (/ = 0) we are lead to 

a = ; q = 1 — p, (6.33) 

or 

- = irTn ' ^ = 2- (6-34) 

a 4(p + l) 

For / 7^ we obtain 

i-P+^{i-py + m + i)p 

a = ; q = . (6.35) 

In case A-2 (0 (r ^ 0) ~ a — br'^, cr > 0) equation fl6.25p becomes 

2/?g ^~,'"~V '?~- + ^(^ + 1) - _ p^,^ ^ g. (6.36) 

ba a 

If / = we must have 

p _ b Ta 

a 2a cr + 2 

or 



q = a + 2, (6.37) 



a = ■ q = (T+l. (6.38) 
For / ^ 

- = 7^/(^ + 1) ; q = 'T, (6.39) 
a 2a 

or 

a = ; g = a + l. (6.40) 

In all these cases the integral of the second order variation of the energy (16. 6p can be 
shown to converge in the limit r ^ 0. 

6.3. Stability of multicomponent scalar solitons. 

The analysis of stability in the scalar case can be extended to the multicomponent scalar 
fields. Now we have integration constants Aj and a degeneration in energy of the SSS 

solutions on the sphere of radius A = \jY^f^i Kf in 3?^. Obviously the variations of the 
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energy vanish for perturbations which remain inside this sphere obtained by modifying the 
constants A, in equation (14.10p in such a way that the "total mean-square charge" A remains 
unchanged. The asymptotic boundary conditions satisfied by the fields 0j(r, A^) (obtained 
from the asymptotic behaviour of the associated one-component scalar field solution (r, A) 
through Eq. fl4.10l) ) are modified by these perturbations and the associated charges (defined 
from Eq. fl4.11l) as Aj) are modified. Consequently, charge conservation condition blocks such 
perturbations and prevents a soliton from evolving spontaneously towards another equal- 
energy configuration in the sphere. 

For general perturbations S(j)i{r) the first-order modifications of the scalar charges take 
the form 



AA,- 



1 



N 



f (Xo)(V50,) - 2 / (Xo) ^ (V0, ■ V6<p,) 



(6.41) 



with X = — Y^iLi 0?- The requirement of charge conservation (AAj = 0) imposes bound- 
ary conditions on the perturbing fields which, as in the one-component case, must vanish 
asymptotically faster than the SSS fields themselves. 

The first-order variation of the energy functional is obtained from the integral of Eq. fl4.5p 
and reads 



^A, 



d-,rV 



f (Xo)<50,V0 



f (^o)V(0) 



(6.42) 



where Eq. (14.12P has been used and a partial integration has been performed. Each term of 
this sum has the form of Eq. (16.41) and vanishes because of the same reasons. Thus the first 
variation of the energy vanishes, which is an extremum condition. The second variation of 
the energy functional takes the form 



Aoe: 



i=l 



Ki=l 



A dr 



and can be written as 



Ase: 



+ 



/(Xo)+2Xo/(Xo)cos2(fi) 
fd6(f), " 



N 

E 

1=1 



/ 4r/(Xo)^ 

i=i 



09 



+ 



dr J 

1 fd5^ 
r'^cos'^{9) \ dip 



(6.43) 



(6.44) 



in terms of the angle fl{r) in the internal space between the vector formed by the radial 
derivatives of the components of the perturbing fields and the direction ^ defined by the 

SSS solution . The second integral in this equation is always positive. If / (Xq) is negative 
the first integral is also positive. Otherwise we have 



/ (Xo) + 2Xo / (Xo) cos2(^]) >/ (Xo) + 2Xo / (Xq) > 0, 



(6.45) 



and, owing to Eq. 



this integral is always positive for admissible many- components 



scalar models with finite-energy SSS solutions. Consequently, all these solutions are statically 
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stable. Moreover, the analysis of the dynamical evolution of small perturbations performed 
for scalar solitons can be straightforwardly generalized to this multicomponent case. Such 
an analysis proves the dynamical stability of these solitons. 

6.4. Static stability of generalized electromagnetic solitons 

Similar analysis of stability can be performed for generalized abelian gauge fields. We 



consider a finite-energy ESS solution I I of the field equations (15.151) {Eo{r),Ho = 0) and 
introduce a small perturbing field {Ei{r), Hi[r)) which does not modify the total electric 
charge of the soliton. The first-order modification of the charge density is obtained by per- 
turbing the first of the field equations fl5.17p . which leads to 

V ■ (T = 0, (6.46) 
where 

^ = + ■ (6.47) 

and the index in the derivatives means that they are calculated for the unperturbed 
solution (note that, owing to the parity invariance, the odd partial derivatives of with 
respect to Y vanish in F = 0). From the integration of (16.461) we see that a must vanish 
asymptotically faster than r~^. Then, the regular perturbing field |£'i(r)| must be damped 
faster than E^ir) itself. This boundary condition is similar to the one introduced for scalar 
models but, owing to the electric charge conservation implicit in the field equations, the 
physical meaning becomes here more transparent. 

Let us now consider the variations of the energy functional under such charge-preserving 
perturbations. The first-order variation is obtained from the integration of Eq. (l5.7p and 
reads 

Aie = -2 / rfgrV ■ \A^B\ + 2 f d?,rA^V ■ a, (6.48) 



where we have introduced the time-like component of the four-vector potential for the so- 
lution (£"0 = — Vv4*^,y4 = 0). This expression vanishes due to the boundary conditions and 
the linearized field equation (16.461) . This is an extremum condition. 

In calculating the second variation of the energy let us expand the first of the field equations 
(I5.17P to the second-order. We are lead to 

V-(a + r/) = 0, (6.49) 
where now the term 



To avoid difficulties related to the gauge determination we work directly with the fields. When the use of the potentials 
becomes necessary in some step of the calculation wc will fix the gauge through appropriate conditions. 
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includes the second-order corrections. Using this equation the second variation of the energy, 
obtained from the integration of Eq.( 15.7p . becomes 



9Xo 
4fV ■ \A^f] 



2 (-^0 ■ Elf 



+ / d^r 



(6.51) 



The last integral in the r.h.s. of this equation vanishes, owing to the boundary conditions. 
The first term is positive if > 0; on the other hand, if < the integrand of this 
term can be written as 



d(p d'^if 



dXo 



'dXl 



(6.52) 



where 6 is the angle between Eq and Ei. By deriving Eq. ( 15.191) with respect to r and taking 
into account the monotonicity of E^lr) we see that the first term in fl6.5ip is always positive. 
Concerning the second term of (16.511) it is positive if < while if > we can write 
its integrand as 

dip 



dXn 



,9V 



{E',cos'i9))]>{H, 



' dip 
'dX. 



(6.53) 



From the arbitrariness of the perturbing fields, the positivity of this term and, finally, the 
positivity of the second variation of the energy requires the condition 



dx - dY^' 



(6.54) 



to be fulfilled in the range of values of X (F = 0) where the ESS solutions are defined. This 
is a necessary and sufficient condition of static stability to be satisfied by the lagrangian 
densities of admissible models supporting finite-energy ESS solutions. This stability criterion 
goes beyond the widely used Derrick's necessary conditions [lOj. 

Let us check, using this criterium, the linear stability of the electrostatic finite-energy 
solutions of the B-I model. From Eqs. (15.51) and (I6.54p we immediately obtain 



dip 
dX'o 



2X, 



OYq^ Stt 



-^(l-/x2Xo)'^'>0, 



(6.55) 



and the stability condition (I6.54p is fulfilled since the ESS field is bounded everywhere 
(Xo < V/x^). 



6.5. Dynamic stability of generalized electromagnetic solitons 

Let us now analyze the dynamical evolution of the small perturbations of the ESS solitons. 
The system of linearized field equations, obtained by expanding (15.171) to first order, is formed 
by Eq. (16.461) aside from a vector equation: 
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V-(T = V-(S-^i) = 

-^-Vx{fl-Hi) = -(S ■ ii) - V X (fl ■ H,) = 0, (6.56) 



where now we have introduced the symmetric tensors 



which will be useful in simplifying the notations in the sequel. The perturbing fields must 
also satisfy the first set of Maxwell equations. Expanding flS.lSp up to the first order we 
obtain 

V-/?i = 0. (6.58) 

Let us look for solutions of these equations which are products of functions of time and 
space variables for both electric and magnetic fields, of the form 

Ei(t,r) = Te(t)-e(f) 

H,it,r) = mt)-h{f). (6.59) 
In this way we are lead to the first-order equations for the time variables 

Th = -f^T„ (6.60) 

where A and fi are separation constants. By deriving these equations we are lead to the 
system 

Te(t)+rTe(t)=0 

n{t)+m{t) = 0, (6.61) 

with r = A ■ /i. The identification A = /i = can be introduced without loss of generality, 
as a consequence of the first-order equations and the positivity of F, which will be established 
below. The eigenvalue F being positive, the solutions (normalized to unity) take the form 

Te{t) = cos{VTt + 6) 

Th{t)=sm{Vft + 5), (6.62) 

where 6 is the same constant phase for both solutions. In this case the eigenfunctions remain 
bounded as time evolves and the soliton is dynamically stable. The field equations for the 
spatial components are 

V-(S-e) = 

V X (O ■ /I) = VfS ■ e, (6.63) 
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where we have used the definitions (16.571) . Note that the first of Eqs. (16.631) is an immediate 
consequence of the second one. Moreover, the first set of Maxwell equations leads to 

V X e = Vfh 

V-h = 0. (6.64) 

We shall now introduce a four-potential A'^ for the perturbing fields, defined in the Hamilton 
gauge {A^ = 0) in such a way that 

H^ = V X A. (6.65) 

This vector potential is determined up to the gradient of an arbitrary time-independent scalar 
field. In terms of this vector potential, the first set of field equations (16.581) are identically 
satisfied while the second set (16.561) becomes 

|(v.(E.i))=0, 

— (S ■ 1) + V X (f2 ■ (V X A)) = 0, (6.66) 

If we write equations (16.651) for the separated functions (I6.59p . by integrating in time the first 
one and using (I6.60p . we obtain the general form of the vector potential for these functions 
as 

A{t, f) = Ta{t)d{f) + V0(f), (6.67) 
where 0(r) is a time-independent function and 

r„(t) = ^ ; a{r) = e{r). (6.68) 

This vector potential, determined in Eq. (l6.67p up to the gradient of a time-independent 
scalar field, becomes univocally fixed by requiring its form to separate in time and space 
variables, taking the form of the first term in the r.h.s of Eq. (l6.67p . In terms of this potential, 
using the definitions (16.570 . the field equations for the spatial part of the perturbation reduce 
to the unique vector equation 

V X (ri ■ (V X a)) = rS ■ a, (6.69) 

which outlines the eigenvalue problem for the linear oscillations in this electromagnetic case. 

As in the scalar case, the standard analysis of this problem can be performed for the 
ESS solitons of admissible generalized electromagnetic field models without any reference 
to the explicit form of the lagrangian density. In this way we shall show that two given 



^''This gauge-fixing condition is allowed by the gauge invariance of the equations for the perturbing fields Il6.56l l- (l6.58t which 
are independent from any gauge choice for the unperturbed fields. 
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eigenfunctions ai(r) and 02 (r) associated to the eigenvalues Fi and respectively, are 
orthogonal and finite-norm with respect to the scalar product 



< ai ■ 02 >= C J d'i'fiai ■ S ■ 02) 



(6.70) 



where C is a normalizing factor. Indeed, multiplying Eg. (16.69 1 ) (defined for a given eigen- 
function ai) by an eigenfunction 02, integrating in space, using the identity 



02 ■ V X (f2 ■ (V X (ii 



0.2 X 



(fl- (V 



X ai 



+ (V X ai) ■ f2- (V X 02), (6.71) 



and taking into account the symmetry of the tensors S and f2, we obtain the equation 



a2 X 



X ai 



+ j dsfiV X ai) ■ ri • (V X 02) 
= '^1 J dzr{ai ■ S ■ 02). 



(6.72) 



Owing to the boundary conditions the first integral in the l.h.s of this equation vanishes 
while the second one converges. By permuting the indices and subtracting we finally obtain 



{Ti - T2) j dsf{ai ■ S ■ 02) = (Fi - r2) < 0.1-0.2 >= 0. 



(6.73) 



We see that the eigenfunctions associated to different eigenvalues are orthogonal with respect 
to the scalar product ( 16.70p . Moreover, if oi = 0,2 we obtain from (16.72^ 



r / dsr{a ■ S ■ a) 



d^fi^V X a) ■ ri ■ (V X a) 



(6.74) 



and both integrals converge. The integrand in the l.h.s. of this equation can be written as 



(a -12 ■ a) 



dX'n 



10] — a 




9F2 



(6.75) 



and is positive for ESS soliton solutions of any admissible electromagnetic model. The inte- 
grand of the r.h.s. of (16.741) takes the form 



(V X a) • ri ■ (V X a 




V X a ) 



- ■ V X a 

r 



- ■ V X a 



+ 



(6.76) 



This term is positive if (and only if) the condition for static stability f l6.54p is fulfilled. 
Under this condition the eigenvalues F are well defined and positive and, consequently, 
the behaviour of any initially bounded perturbation remains bounded as time evolves. We 
conclude that the statically stable ESS solitons of admissible electromagnetic models which 
satisfy fl6.54p are also dynamically stable. Moreover, the spectrum of eigenvalues is discrete 
and the eigenfunctions generate the functional space of the physical vector potentials, which 
can be written as 

A{t,f) = ^C„sin(v^t)a„(f), (6.77) 
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and are in a one-to-one correspondence with the physical perturbed states of the sohton. 
Indeed, any charge-preserving perturbation of the sohton field is described, in the Hamilton 
gauge, by vector potentials which can be obtained from one (and only one) of the form f l6.77p 
by the addition of gradients of time-independent scalar functions. 

The analysis of the spatial structure of the eigenfunctions and physical perturbations can 
now be performed by separating in radial and angular parts the components of the vector 
potentials an{r) in the natural basis of the polar coordinate system. This procedure, which 
is standard in spherically-symmetric physical problems [51j, will determine the asymptotic 
and central-behaviour of the perturbing fields, as in the scalar case already considered. We 
shall leave this study for future developments. 



6.6. Stability of generalized non-ahelian gauge solitons. 



Owing to the essential self-interactions involving the field potentials, the treatment of 
the static and dynamic stability of the solitons for generalized non-abelian gauge models is 
more involved than in the abelian case. A detailed study of the stability of some extended 
static finite-energy solutions for the standard Yang-Mills model has been performed in Ref. 
^47j . where similar difficulties arise. Most methods of that work can be generalized to the 
present situation and we shall follow this way in analyzing the stability behaviour of the 
ESS solutions fCTD . 

Consider a finite-energy ESS solution of the field equations (15.331) of the form (15.401) and 
introduce small regular perturbing fields through the definitions 

Ea (r) = -VA°(r) , #a = , 1, = 

5Ea{f) , 5Ha{^ , M°(f) , M„(f). (6.78) 
To first order these fields are related through 

be 

6Ha{f^=V x5Aa, (6.79) 

and they are assumed to leave invariant the color charges Qa associated to the unperturbed 
solution. To first-order the modifications of the charge densities are obtained from the per- 
turbation of the time-components of the field equations (I5.33P (the generalized Gauss laws) 
and read 

V ■ = -^7 E Cabc^SA, ■ E,, (6.80) 

be 

where 

dip 



(7a 



(9X, 

and the modifications of the total charges read 
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Air J J 9Xo 



(6.82) 



The first-order perturbations of tlie vector equations, given by tlie spatial components of 
Eqs.f l5.33|) (tlie generalized Ampere laws) read 



Ot 



where 



dip 
dXc 



(6.83) 



(6.84) 



The r.h.s. in Eq.f l6.80|) is the color-charge density carried by the perturbations and its spatial 
integral must vanish, according to our initial assumptions. This requirement restricts the 
asymptotic behaviour of <?„ and, as for the other field models already considered, leads to 
boundary conditions to be satisfied by the perturbing fields. Moreover, similarly to the multi- 
component scalar case, the color charges of the unperturbed solution (Qa) fix a direction in 
the color-charge space (called in Ref.[17j "electromagnetic" direction, while the orthogonal 
directions are termed "charged"). Owing to the first-integral equation (15.381) . the potentials 
and the fields Ea of the ESS solutions lie in this direction. For the perturbing fields 
to remain purely electromagnetic the associated charge densities that are induced by them 
must vanish and, owing to Eq. (l6.80p . 6Ab must also lie in this direction. In what follows, we 
shall prove the stability of the finite-energy ESS solutions against this kind of non-charged 
perturbations. 

Now let us analyze the variations of the energy functional. The first variation is obtained 
by perturbing the spatial integral of fl5.35p (this is a gauge-invariant quantity, as well as its 
variations) around the ESS solutions, which reads 



Aie = -2 / dsfJ2 V ■ \A%] + 2 / dsfJ2 V ■ a,. 



(6.85) 



The first integral in this expression vanishes, owing to the boundary conditions. Using 
Eq. (16.801) the variation becomes 



Aie = -2g [ d.fJ^Catc^A'idA, ■ E,. 

abc 



(6.86) 



As expected this expression vanishes, owing to the parallelism of A° and E^. in the color 
space (see Eqs.f l5.40|) and fl5.43p ) and the antisymmetry of the structure constants. 

In obtaining the second variation of the energy functional we follow the same steps as in 
the abelian case. First we expand the generalized Gauss law to the second order. After the 
cancellation of the first-order terms we are lead to 

9V 



V ■ Qa 

where now 



-g X! ^<^bc 

be 



dip 
dXo 



5A, ■ 5E, + 2^ [Y^E,-5Ep] 5A, ■ E, 



(6.87) 
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Co-, 



9V 



+ 2 



p 

3, 



5Xo \ p 



(6.88) 



By expanding the integral of fl5.35p up to second order and using Eqs. (16.801) and (I6.87P 
the second variation of the energy becomes 



A2e = y d^r 



+ 



dip 



(6.89) 



- 2 / d'fV ■ (AS ■ cD'^) -2g f d,r^ Yl CabcA'jA, ■ 5E, 



ahc 



Once again, the divergence term in this expression vanishes owing to the boundary con- 
ditions. The integrand of the last term is the scalar product in color space between the 
potential of the unperturbed field and the first component of the second-order perturba- 
tion of the color-charge density in the r.h.s of Eq. fl6.87p . For electromagnetic perturbations 
this component must satisfy the condition 



dip 



Y,Cai,cAl5A,.5E, = Q, 



(6.90) 



abc 



(note that the remaining component in the r.h.s. of Eg. (16.870 lies already in the electro- 
magnetic direction). Consequently, the last term in Eg. (16.890 must vanish and the second 
variation of the energy takes a form similar to that of the abelian case. We can now deter- 
mine the conditions for stability of the finite-energy ESS solutions in this non-abelian case 
through a similar argumentation. As easily seen stability requires the lagrangian-density 
function ip{X^ Y) to satisfy the condition 



dx - dY^' 



(6.91) 



in all the range of values of the gauge invariants (X, Y = 0) defined by the solution. This 
condition is formally the same as in the abelian case and is also necessary and sufficient 
for the stability of the solitons against electromagnetic perturbations. Obviously, it is a 
gauge-invariant criterion. 

The analysis of the dynamical stability of non-abelian solitons should now be performed 
starting with Eqs. ( 16.801) and ( 16.830 and following similar steps as in the abelian case. But 
the presence of the antisymmetric structure constants and the symplectic character of the 
eigenvalue problem require new qualitative procedures and longer calculations which would 
lengthen excessively the contents of this work. This issue will be approached elsewhere. 
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7. Conclusions and perspectives 



In this work we have solved the problem of characterizing a large class of physically 
consistent relativistic lagrangian field theories in three-space dimensions, supporting static 
spherically symmetric non-topological soliton solutions. The fields concerned were one and 
many-components scalar fields (whose lagrangian densities depend on the kinetic term alone) 
and generalized gauge fields of compact semi-simple Lie groups. This characterization is 
exhaustive and leads to the classification of such models into six types, according to the 
central and asymptotic behaviours of the soliton fields. We have performed a broad analysis 
of the linear stability of the solutions, obtaining necessary and sufficient stability conditions 
which go beyond the usual Derrick criterion. We also have carried out a general spectral 
analysis of the linear perturbations around the soliton solutions, confirming their dynamical 
stability and setting grounds for their quantum extensions. 

All these results allow the explicit determination of a large number of examples of such 
a class of lagrangians, providing a wide panoply of tools for the analysis of diverse physical 
problems, as those mentioned in the introduction and others. Among these problems let us 
outline three of particular interest, which we are addressing from the methods developed 
here. 

1) As already discussed in section 3.1., the photon-photon interaction mediated by the 
QED vacuum can be classically described in terms of effective lagrangians which are poly- 
nomial expressions in the gauge invariants X and Y and can be obtained in a perturbative 
procedure [21], where the lowest order is the well-known Euler-Heisenberg lagrangian ^32j 



{rj and ^ being positive constants). The sequence of these lagrangians exhibits finite-energy 
ESS solutions and suggests an interpretation in terms of the screening effects of the vacuum 
on the field of point-like charges. Unfortunately, as already mentioned, the perturbative 
expansion involved in this procedure is a low-energy (or a low-intensity field) approximation 
and is not accurate to describe the strong fields present near the center of the ESS solutions. 
It is thus necessary to explore this issue with other effective lagrangians, obtained from 
the perturbative renormalization of the self-energy of point-like fields. The analysis of this 
problem is in progress ^T], |52j . 

2) Scalar field models as that of the example treated in section 3.4., which belong to 
the case B-3, exhibiting short-ranged SSS solutions (solitons or not), can be extended to 
generalized gauge field models supporting similar ESS solutions. This behaviour, which is 
related to the form of the lagrangian density around the vacuum, may arise in effective 
actions for some fundamental forces. If we assume the effective dynamics of the non-abelian 
gauge fields in electroweak interactions to be described by this kind of lagrangians, the short 
range of these forces could be explained in terms of the non-linear self-couplings among 
these fields, coming from the integration of some higher energy degrees of freedom of a more 
fundamental theory. In this case the appeal to any symmetry breaking or Higgs mechanism 
should become superfluous. In our sense, this alternative deserves to be thoroughly explored. 

3) A new approach to the phenomenological description of the hadronic structure can be 
envisaged, using the results of section 6 on the spectral analysis of the excitations of the 




(7.1) 
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(multi-) scalar or generalized gauge-invariant solitons. In the phenomenological descriptions 
based on the Skyrme model the hadron arises as a topological soliton of a non-linear field 
theory. Other models, which are believed to give an effective low-energy approach to the 
non-perturbative regime of QCD (as the Friedberg-Lee model [^,[8] and related theories), 
describe hadrons as confined states of quarks in non-topological solitonic bags of non-linear 
phenomenological fields. As an alternative to these viewpoints, we are considering a gener- 
alized gauge-invariant lagrangian model for gauge fields ("gluons") coupled to a quark-like 
fermionic sector and implementing properly chosen symmetries. Such a phenomenological 
model may be interpreted as an effective lagrangian for QCD or, alternatively, as a field- 
theoretical low-energy limit of string theory. The classical generalized gauge-field lagrangian 
can exhibit soliton solutions in absence of other fields. If such solutions are minima of the 
functional of energy of the full action, their small perturbations will involve fermionic and 
bosonic modes. The quantization of these modes leads to "quasi- quarks" and "quasi-gluons" 
as quantum excitations of the soliton field. This quantum extension becomes a model for the 
hadron containing these particles. In this picture the confinement would be a consequence 
of the fact that quarks are quasi-particles associated to these quantum excitations and (as 
the phonons in a solid) they cannot exist outside the hadron. 

Another domain where the results of this work could be useful concerns the search for 
self-gravitating (scalar and gauge) field configurations in General Relativity [23j. The classi- 
fication of the lagrangian field theories considered here, supporting non-topological solitons 
in fiat space, can be extended to the static spherically symmetric solutions of the Einstein 
equations resulting from the coupling of these fields to gravitation. Indeed, we have verified 
that these equations have first-integrals which have the same form as (or can be closely re- 
lated to) the ones obtained from the corresponding field theories in flat space (of the generic 
form of Eqs. fi2.4p or fl5.19p ). This result opens the possibility of generalizing to the gravita- 
tional case many of the methods and results obtained here [53j- We will continue to address 
this topic in future work. 

It would be also interesting to study the soliton solutions of generalized non-abelian gauge 
field theories with other ansatzes than the ESS one. In fact, such a kind of solutions have 
been already found for the SU(2) non-abelian B-I theory within the "monopole ansatz" [TJ]. 
This issue should be a theme for a future investigation. 

Let us conclude with some comments concerning an important question which has not 
been considered here. It refers to the analysis of propagation of wave-like solutions of these 
models. As can be easily seen, all these theories exhibit plane wave solutions propagating 
with the speed of light. But, owing to the non-linear self-coupling, they also support other 
radiative solutions propagating with more complex dispersion relations. In most cases such 
waves evolve towards spatially-singular configurations. Roughly speaking, the wave fronts 
travelling with velocities which are dependent on the values of the fields at every point tend 
to cumulate, generating discontinuities and shocks after a critical time. Regularly evolving 
wave solutions of a system of field equations are called exceptional. If all the wave solu- 
tions of a given system are exceptional, the system is called completely exceptional [54]. A 
detailed analysis of the problem of wave propagation for generalized electromagnetic field 
models was performed by G. Boillat [55], who established the complete exceptionality of 
the Born-Infeld electrodynamics. Moreover, B-I is the only admissible generalized electro- 
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magnetic field theory (with asymptotically coulombian elementary solutions) exhibiting this 
property. Nevertheless, the Boillat analysis considers only models which satisfy the condi- 
tion ^{X = 0,Y = 0) = 1 (case B-2) and, consequently, excludes the models belonging 
to B-1 and B-3 cases. It would be interesting to extend the Boillat analysis for these cases. 
However, when one considers the extensions of B-I electrodynamics to the non-abelian case 
or in the Kaluza-Klein context, this exceptionality character is lost [56] . 
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